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EXTENSIONS TO MISSING PLOT TECHNIQUES 


H. R. THompson 


Applied Mathematics Labo;atory, Department of Scientific and Industrial Research 
Wellington, New Zealand 


When several plots in an experiment have missing observations it 
often Lecomes very laborious to apply the simple method of giving 
preliminary estimates to all except one, calculating it by the formula 
for a single missing plot and proceeding iteratively through two or more 
cycles. The technique set out in this note permits the simultaneous 
estimaticn of missing observations in a randomised blocks experiment 
with r blocks and ¢ treatments, for the special case when there is one 
missing observaiicn in eac’: of n treatments (n < #), distributed over 
p blocks such that n; are in the ith block. (The entirely analogous case 
of one missing observation in each of % blocks (n < r) is obtained by 
interchanging r and ¢ in the succeeding formulae.) More complicated 
eases, e.g., some treatments having more than one missing observation, 

‘ean be solved iteratively, either by a combination of the single plot 
technique and the present technique or by applying the present tech- 
nique twice, depending on the distribution of the missing observations. 
The method is likely to prove useful when the data from part of one or 
more blocks is missing, due to reasons beyond the experimenter’s 
control. 

Let the missing observations in the ith block (i = 1, --- , p) be 
denoted by z,;(j = 1, --- , n,), the totals of existing plots in the treat- 
ments corresponding to these observations by 7%; , the total of existing 
plots in the ith block by B/ and the grand total of existing plots by G’. 
The simultaneous estimation by least squares of the 2z,; leads to the 
equations 


fi-1 
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241 


ar 
‘ 
{ 
} 


242 BIOMETRICS, SEPTEMBER 1956 


which are, in matrix notation, 
=C, 


where x is the column matrix (z;;), C is the column matrix 
(rB; + — G’), and 


Ka, A, K., 

K,, K, A, 


in which A; isan; X n; matrix with elements a; = (r — 1) (¢ — 1) in 
the main diagonal and b; = —(r — 1) elsewhere, 
K,; isan, X n; matrix with all elements equal to k;; (k;; = 1 here). 
The solution of the normal equations may be shown to be 
z=J''C, 
where 


Kip 


say, 


L a Ki, A 
in which the A{ , Kj; are also A- and K-type matrices respectively, with 
elements a{ , bj , ki; . Given below are the solutions for p = 1, 2, 3. 


p=l. a —n) 
= 1/(r — 1)t(t — n) 
p=2. af = — — m+ — m) — — 
by = [¢ — 1)°(t — + 
az = [(r — 1)°(t — n,)(t — m + 1) — mlm — 1)J/D 
= — 1)°(t — m) + 
kin = kh, = — 1)t/D 


ll 


where 


D = — — — — m2) — muna]. 
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p= 3. Writing C,=(¢ I(t—n), 

ay = (Ci +r — — nama) 

— (rm — + — 2nsns)]/D, , 
bi = — 1)(C.Cs — nana) + + — 2nsns)]/D, , 
= — IMC, — , 

kis = — — , 
where = (r — 1)t[C,(CsCs — mate) — miltaCs + — 
and similarly for the other two sets by interchanging the suffices. 


EXAMPLE 


In an experiment on virus prevention in swedes (for the data of 
which the author is indebted to Mr. P. R. Fry, Plant Diseases Divisign, 
D.S.1.R.) with 14 treatments and 6 blocks, there were 9 missing observa- 
tions, one in each of 9 treatments, distributed over 3 blocks as indicated 
in Table 1, in which the various treatment and block totals are also 


TABLE 1 
Missing Plot Calculations 
Treatment By G’ 6B,’ + 14T,;' 
10 298.9 666.8 3697.5 4487.9 62.1 
2 169.7 566.8 ? 2079.1 39.2 
ll 215.1 a 4 2714.7 48.3 
13 211.9 2669.9 47.6 
14 280.4 3628.9 61.3 
1 392.1 549.8 . 5090.7 87.4 
3 221.4 2700.9 53.2 
4 239.7 2957.1 56.9 
12 250.7 3111.1 59.1 


given. The original data gave the percentage infection in each plot, 
and these percentages have been transformed to equivalent angles. 
Substituting r = 6, ¢ = 14, n, = 1, n. = n; = 4 in the equations for 
p = 3, we find 


a 
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| 
2819110 44345 — 895 
—805 44406 
—805 4133 
—805 4133 
—805 4133 
—805 —1120 
—1120 
—805 —1120 
(-805 —1120 


The calculations can be quite easily systematized to obtain the z,; , 
and these are given in the last column of Table 1. It is interesting to 
note that if a preliminary estimate of the first missing value is made 
and the other 8 calculated using p = 2, the values obtained after one 
cycle are all within 0.1 of their “correct’’ estimated values, whereas if 
the solution for p = 1 is used and preliminary estimates made of the 
first five, two cycles are required to obtain the correct values and the 


4133 


—1120 


—1120 


—1120 


—1120 
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—805 —805 
4133 4133 
4133 4133 

44406 4133 
4133 44406 
—1120 »—1120 
—1120 -—1120 
—1120 -—1120 
—1120 —1120 


calculations are more laborious. 


4133 


4133 


4133 


4133 
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| 
—805 —805 -805 -805 —805 
4133 —1120 -1120 -1120 —1120 
44406 -1120 -1120 -1120 —-1120 
4133 —1120 -1120 -1120 —1120 
-1120 -1120 -1120 —1120 
44406 4133 4133 4133 
4133 44406 4133. 4133 
| 
i 
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THE ANALYSIS OF A 3 X 6 EXPERIMENT ARRANGED IN 
A QUASI-LATIN SQUARE 


G. E. Hopnetr 
Rothamsted Experimental Station, Harpenden, Herts., England 


Introduction 

A factorial experiment on Star Grass (cynodon sp.) to investigate the 
effect of three frequencies of cutting and six levels of nitrogen was 
commenced in 1953, at the Grasslands Agricultural Research Station, 
Marandellas, S. Rhodesia. The design was derived from that given by 


TABLE 1* 
Design, Field Plan and Yields 1953/54 
[Dry matter: lb. per plot (1/200 acre)] 


Ce J; Ji Ji J: J: J; 


J; Js J; J; Ji ji Total 
I; h 200 221 120 lll 010 001 
_3.5 28.6 19.3 29.9 17.8 48.6 147.7 
1; i 221 210 111 100 001 020 
31.9 9.7 30.5 3.5 38.4 34.0 148.0 
I; qT; 020 011 210 201 100 121 
31.0 51.4 9.1 22.9 2.2 30.3 146.9 
ih I; 110 101 000 021 220 211 
10.5 23.2: 13.0 55.4 16.2 29.7 148:0 
I, ds 101 120 021 010 211 200 
24.3 21.8 57.1 17.4 24.5 4.3 149.4 
nh od O11 000 201 220 121 110 
42.3 4.4 23.5 16.0 41.6 16.7 144.5 ~ 


Total | 143.5 139.1 152.5 145.1 140.7 163.6 884.5 


200 = asboco, 221 = asbyti, etc. 


*The Z and J components of the interactions AB and ABC are partially confounded respectively 
with the rows and columns of the square, the J's and J's above being the diagonal sets of combinations 
of a and b as defined in Table 2 (iii). 
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Yates [1] fora 3 X 3 X 2 experiment (factors A, B, C) arranged in a 
6 X 6 quasi-Latin square, by formally identifying the six levels of 
nitrogen with the six combinations of the pseudo-factors B and C as in 
Table 2. The design, layout and yields in the first. season are shown in 
Table 1. 

(Treatments: Frequency of cutting (a): 2, 3, 6 cuts at regular intervals 
over a period of 18 weeks during the rainy season. Nitrogen (b and c): 
0, 500, 1000, 1500, 2000, 2500 lb per acre sulphate of ammonia (21% N), 
each level being applied in sixths at intervals of 3 weeks. Basal dressing: 
400 lb per acre superphosphate with the first dressing of nitrogen.) 

~ As there was considerable interest in the interaction between the two 
factors, the two-way table of results was required free from row and 
column effects, together with the appropriate standard errors. In order 
to make the adjustments, it was necessary to evaluate’ the formal 
three-factor interaction, since it was a component of the interaction 
between nitrogen and the frequency of cutting. In the course of this 
analysis, a general method was developed for determining the standard 
errors of comparisons between adjusted yields. This note puts on record 
the methods used in the analysis of this unusual design. The experiment 
is being continued and the agricultural implications of the results will 
be discussed elsewhere [2]. 


Analysis 


The computations were set out asin Table 2. The sums of squares for 
rows and columns (ignoring treatments) and for the main effects of the 
two factors were computed in the usual way from the marginal totals 
of Tables 1 and 2 (i) respectively. The table of sums (c, + c)) was the 
starting point for computing the effects of A, B and AB while the inter- 
actions of these effects with C were computed from the table of differ- 
ences (c, — ¢ ). For the main effect of C, the sum of squares was 
obtained from the total of the latter table as (383.7)’/36. As a check, 
the total of the sums of squares for the components B, C and BC was 
compared with the sum of squares for nitrogen already obtained. The 
linear effect of nitrogen N’ was also evaluated, as the nitrogen totals 
indicate a very marked linear response. 

For the interaction AB, the estimate of the components and the 
corresponding sums of squares were computed as in Table 2 (iii). The 
ordinary J and J totals of the (c, + co) table were adjusted by means of 
the row and column totals (R, --- R, , C, --- Cs) in order to eliminate 
the effects of rows and columns respectively. These expressions for 
the adjusted totals (J’, J’) may be derived by the standard method of 
fitting constants (Yates [3], Nair [4], Kempthorne [5]). Substitution 


i 
2 


QUASI-LATIN SQUARE 


of the constants gives 
21; = 244,+ 612+ 6%; 
213 61, + 247, + 
21; = 6i,+ 67%. + 24%; 


which also verifies that the I’ are free from row constants. Although 
the J’ are correlated, differences between them are proportional to the 
differences between the constants. For example: 


— 12) = 18(4, — 


so that 
~ 
1 2 


Thus the estimates of the J component of AB (2 d.f.) are given by the 
deviations of the I’ from their mean divided by 9 (computed as A2/’/18). 
In the absence of confounding the corresponding estimates would be 
AI/12, so that the partial confounding has reduced the ‘effective’ 
replication from 12 to 9 i.e. by a factar of 3/4 (the relative information). 
The normal divisor for the sum of squares (here 12) is also modified 


TABLE 2 
Computations 
(i) Treatment totals 
Nitrogen 
Frequency No m ns 
of Total 
cutting Co Co Co 
bo bh be bo bs 
a 17.4 35.2 65.0 87.0 93.7 | 112.5 | 410.8 
a 5.7 27.2 41.1 47.5 60.4 71.9 | 253.8 
a, 7.8 18.8 32.2 46.4 54.2 60.5 | 219.9 
Total 30.9 81.2 | 138.3 | 180.9 | 208.3 | 244.9 | 884.5 
Linear com- 
ponent N’ —5 -3 -1 +1 +3 +5 
S.S. for N’ = 70 (-—5 X 30.9 —3 X 81.2... +5 X 244.9)? = 5313.66 
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G — Ce 
Total | be b, b, | Total 
410.8} 69.6 | 58.5 | 47.5 | 175.6 
253.8 | 41.8] 33.2] 30.8 | 105.8 
219.9 | 38.6] 35.4] 28.3 | 102.3 
884.5 | 150.0 | 127.1 | 106.6 | 383.7 
Al = AI/12 
—10.1 —0.84 
+ 8.9 +0.74 
+ 1.3 +0:11 
j = 4J/12 
— 4.4 —0.37 
—20.0 —1.67 
+24.5 +2.04 
Adjusted component A2I’ = A2/'/18 
21, = 21, + Ri + Rs = 868.7 —18.8 —1.05 
21, = 21, + Ri + Re = 899.6 +15.1 +0.84 
21,’ = 21, + R, + Ry = 888.2 + 3.7 +0.21 
Total 2653 .5 
Mean 884.5 
Adjusted component A2J' = A2J'/18 
= 2W,4+C,+ C, = 869.4 —15.1 -—0.84 
= 23+ Ci + Ce = 865.7 —18.8 —1.04 
= 23,+C, + C, = 918.4 +33.9 +1.88 
Total 2653.5 
Mean 884.5 


8.8. for AB = Po (18.8 + 15.1% + 3.7% + 15.1° + 18.8 + 33.0%) = 64.61 


(ii) Two-way tables 
| 
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(iv) From a — Ce 
Component 4H h = AH/12 
HA, = 131.1 +3.2 +0.27 
Hy = 124.7 —3.2 —0. 
H,; = 127.9 0.0 0.00 
Total 383.7 
Mean 127.9 
Component 4K k = AK/12 
K, = 135.8 +7.9 +0.66 
= 128.6 +0.7 +0.06 
K,; = 119.3 —8.6 —0.72 
Total 383.7 
Mean 127.9 
Component A2H’ h! = A2H’/6 
2H,’ = 2H, + Ri — R: — Re + Re = 259.0 +3.2 +0.53 
2H,’ = 2H: +R, + Rs — — BR: = 246.9 —-8.9 —1.48 
2H,’ = 2H, + Ri — Ri + Rs — Re = 261.5 +5.7 +0.95 
Total 767.4 
Mean 255.8 
Component A2K’ k’ = A2K’/6 
2K,’ = + C; Cs Cs = 258.9 +3.1 +0.52 
2K,’ = 2K; C; C; Cs Cs = 260.4 +4.6 +0.77 
2K,’ = C; —C,—C.+ Cy = 248.1 -—7.7 —1.28 
Total 767.4 
Mean 255.8 


8.8. for ABC = = 3.2 + 8:9 + 5.79 + 3.12 + 4.6 + 7.7%) = 17.67 


by this factor. Thus the sum of squares for AB (J) is 


1 
= 36 S(A2I’) 


which is the form most suitable for computations. The J component is 
treated similarly. 


Section (iv) of Table 2 shows the corresponding computations for the 
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formal three factor interaction ABC, where H and K are the J and J 
totals of the (c, — c ) table. In this case, for example 


hy — = 3 (Hi HD 
whereas in the absence of confounding 
h, hy = — H,) 


so that the relative information on ABC is (3)/6 = }. By the rule 
used above, the sum of squares for ABC (J)(2 d.f.) is given by 


1 __ 1 
x 12) = S(A2H’) 
TABLE 3 
Analysis of Variance 
Source of variation D4. Sum of squares Mean square 
Rows 5 2.24 
Columns 5 70.45 
Frequency of cutting 2 1728.92 864.46 
Nitrogen 5 5395.80 1079.16 
linear (N’) 1 5313.66 5313.66 
remainder 4 82.14 
Interaction 10 367.20 
AB 4 64.61 
AC 2 284.92 
|ABC 4 17.67 
Error 8 202.27 25.28 
Total 35 7766.88 


S.E. per plot, s = +5.03 or 20.5% 


The complete analysis of variance is shown in Table 3. Both the 
main effects were clearly significant, the sum of squares for nitrogen 
being almost entirely due to the linear component... Although the inter- 
action mean square (10 d.f.) did not reach significance, the magnitude 
of the AC component indicates that there was some effect. From the 
table of adjusted mean yields (Table 4), the response to nitrogen is 
readily seen to have been reduced as the frequency of cutting was 
increased. The nature of this interaction is examined in more detail 
later. 
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TABLE 4 
Treatment Means Adjusted for Confounding 
(Dry matter: lb. per acre) 


Nitrogen 
Frequency No Ne Ns 
of Mean 
cutting Co Co Co CQ CQ Cy (+1. 45) 
bo | b; be bo b be 
Qo 7.9(a) | 16.7(b) | 34.2(c) | 43.0(j) | 49.2(k) | 54.4() | 34.2 
a 4.1(d) | 13.5(e) | 19.4(f) | 24.0(m) | 29.5(n) | 36.4(0) | 21.2 
3.5(g) | 10.4(h)| 15.6(2) | 23.5(p) | 25.4(q) | 31.7(r) | 18.3 
Mean (+2.05) | 5.2 13.6 23.1 30.1 34.7 40.8 24.6 


The adjusted treatment means (Table 4) were obtained using the 
following formula: 


= Jae +t’ + 
+ (h’ + k’ —h — k) 


where the required adjustments are determined by the J and J sets to 
which the treatment combination belongs; the plus sign is used when the 
treatment.is at the c, level and the minus sign when it is at the cy level. 
Consider, for example, the adjustment of the mean for adoboco . This 
treatment belongs to the sets J, and J, and is at level c,. Thus from 
Table 2 

17.4 


— (0.53 + 0.52 — 0.27 — 0.66) = 


Computation of standard errors of comparisons between adjusted treatment 
means 


Because of the loss of information on some of the components of the 
interaction due to the confounding, the standard errors of comparisons 
within the body of Table 4 will be increased by factors which depend 
on the particular components of interaction involved. The standard 
error of any desired comparison may be computed quite simply using 
the scheme shown in Table 5 which will be described step by step. 
The underlying theory is given in the appendix. 
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The 17 orthogonal degrees of freedom representing main effects and 
interactions of the factors A, B, C are written down in terms of the 
treatment combinations, which are labelled a, b, c --- g, r in Table 4. 
For example, the linear effect of B i.e. B’ is 


the effect of the other pseudo factor C is 
G+b+ --- +9) --- +9 


and so on for the quadratic components and the interactions. The 
relative information (R.I.), given by Yates [1], on each of these contrasts 
is also noted. The sums of squares (S.S.) of the numbers in each row 
are then tabulated (Table 5). 

If the variance of a simple difference between two mean yields, e.g. 
(q — a) is required, then the numbers in the column headed ‘a’ are 
subtracted line by line from those in the ‘g’ column. The results are 
shown in the column headed (q — a). The factor by which the normal 
variance of a treatment mean (here 25.28/2, since there are two repli- 
cates) must be multiplied is then given by the sum of the quantities 


(q — a)” 
(S8.) R.1.) 


which are computed line by line. For the difference (q — a) this factor is 


2° (-3)' _1)° 


so that var (¢ — a) = (34/9)s’/2 = 47.75. As a check the factor becomes 
2 for a comparison in which the relative information is unity for all the 
orthogonal contrasts involved. 

The variances of the different types of comparison between pairs of 
treatment combinations are given in Table 6 together with examples and 
the numbers of each type. The average variance of the comparisons can 
then be computed if required and hence the standard error which can be. 
assigned to the treatment means. In experiments where there is little 
loss of information the average standard error will be satisfactory for 
most purposes. 
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TABLE 6 
Type of comparison Example | No Variance # 

Factor levels differing 
A. d-a 18 
B b-a 18 Ae 
A and B h-d 36 
Cc n—e 9 
‘Cand A r—c 18 
Cand B k-a 18 43 
CandAandB | q-a 36 az 
Total or mean 153 13 X 21 
17X9 


As the response to nitrogen was substantially linear, it was decided 
to examine the interactions of the other factor with this’ component 
ie. A’N’ and A”N’. The values of these components cofaputed from 
Table 4 were —2.07 + 0.60 and +0.50 + 0.35 respectively in units of 
a single plot. The interaction A’N’ is thus significant. 

The computation of the standard errors of these components provides 
an example of how the above method can be readily extended to cover 
all types of linear comparisons between the adjusted treatment means. 

Writing A’N’ in the form 


the value of each term may easily be computed line by line (as for 
(q — a) above) and afterwards combined. The results are shown in 
the last column of Table 5. Applying the above rule __ 


_ 452 


2 
6 


If the fraction 4/3 is replaced by 1, the variance becomes 70 X s’ as it 
should be if there was no confounding. Similarly var (A’’N’) is found 
to be 4528/2. 

This is actually a rather simple case since only two orthogonal 
contrasts are involved. It is easy to verify that 


4 
| 
4 
| 


QUASI-LATIN SQUARE 255 


N’ = 2B’ + 3C 
so that 
A'N’ = 2A’B’ + 3A'C 
and 
A"'N’ = 2A"B! + 3A’’C 
Hence 


var (A’N’) = 4 Var (A’B’) + 9 Var (A’C) 


4 
= as before. 


This example illustrates the principle underlying the above procedure: 
the required comparison is expressed linearly in terms of the 17 ortho- 
gonal contrasts, whose variances are readily available. Since the 
contrasts are orthogonal, the required variance is obtained as the sum 
of the-variances of the separate terms. 

The method can also be readily applied in the computation of stand- 
ard errors of comparisons in all types of split-plot experiments. The 
error mean square divided by the relative information in the confounded 
experiment is merely replaced by the appropriate error mean square 
(on a sub-unit basis) from the split-plot analysis. 


Acknowledgements 


* My thanks are due to Miss 8. Cohen for her assistance in the analysis 
of this experiment and to Mr. H. D. L. Corby, Grasslands Agricultural 


Research Station, Marandellas, S. Rhodesia for permission to use the 
data in this paper. 


This work was carried out :nder a Colonial Development and Welfare 
Scheme. 


Summary 
The analysis of a 3 X 6 factorial experiment : anged ina 6 X 6 


quasi-Latin square is described and details are given of a general method 


for computing the standard errors of comparisons between treatment 
means adjusted for confounding. The application of this method to the 


computation of standard errors of treatment comparisons in split-plot 
designs is indicated. 
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APPENDIX 


(a) The problem is to express any given comparison as a linear 
function of the orthogonal contrasts, i.e. we require the values of the 
coefficients A, B --- Q, R, of these contrasts. 

Using matrix notation and referring to Table 5 let 


| A) G =.the 18 X 18 array of 

b B coefficients defining 

i the orthogonal contrasts 
x= u= 

q Q 


i.e. z is the vector of treatment combinations and u is the vector of 
coefficients of the orthogonal contrasts. 

The set of orthogonal contrasts: A’, A’, B’ --- (including the total) 
is given by 


Gx 


Now let ¢ be a column vector of coefficients defining that comparison 
whose error is required, so that the actual comparison is given by 


t’x (2) 


In order to express this contrast as a linear function of the A’, A”, B’ 
etc., each of these is multiplied by a coefficient A, B, C --- Q, R and 
identified with 

Thus 


u’Gx = t’x (3) 
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and so, comparing coefficients of x 
uwG=t’ (4) 


* Since the rows of G are orthogonal, GG’ is a diagonal matrix, and 
hence the inverse of GG’ is also a diagonal matrix whose diagonal 
elements are the reciprocals of the corresponding elements of GQ’. 
The diagonal elements of GG’ are given by the sums of squares of the 
elements in the corresponding rows of G; they are given in the column 
headed §S.S. in Table 5 and may be denoted by (SS), , (SS)e -++ (SS)p. 

Postmultiplying (4) by G’ gives 


u’GG’ = t’G’ = d, say 
where 
d = (d, ds dp) 


Hence 


uw’ - ( d, dz dp ) (5) 
(SS), (SS)s  (SS)r 
which is the required solution. 
(b) We now require the variance of the comparison given by ?/’z: 
var (t’x) = var (u’Grx) from (3) 
= wu’ var (Gx) u 
= B’ var (A’) + C’ var (A”) + --- + R’? var(A”B’C) (6) 


since var (Gx) is a diagonal matrix, the contrasts Gx being orthogonal. 
Now the variance of any contrast X(X = A’, A’’, B’ etc.) is given by 


(7) 
where r is the replication. 
Thus 
(S.S.)x 8 
__dy_ 
(8) 
which is the form used in the text. 
(c) Examples 


(1) Suppose var (q — a) is required. Then 
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= (—-100000000000000010) 

d = t/G’ = (0201 —32 —11 —302 —1 —11 —1) 
i.e, the elements of d are obtained by subtracting row by row ‘column @’ 
of Table 5 from ‘column q’; these elements are given in the column 
(q — a) of Table 5. Substitution in (8) immediately gives var (q — a) = 
(34/9)s/2 as shown in the text. 

(2) For the comparison A’N’ 
t = (531000 —5 —3 —1 —1 —3 —500013 5) 
d= = (000000 16000360000000) 
Then 


uw = (00000044 00080000000) 


using (5) and so A’N’ = 2A’B’ + 3A’C as quoted. Var (A’N’) then 
follows using either (8) or (6) and (7) as shown in the text. 
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A NOTE ON THE 4" SERIES OF FACTORIAL EXPERIMENTS 


P. J. CLARINGBOLD 


Department of Veterinary Physiology, University of Sydney, 
N. 8. W., Australia 


In the initial stages of an experimental investigation one of the 2” 
series of factorial experiments is ideal for quick evaluation of important 
factors and interactions. In subsequent and more detailed studies of 
the response surface, designs involving more than two levels of factors 
are often required. The principles involved in fractional replication 
and confounding of factorial experiments are well established [Yates, 
1937; Fisher, 1942, 1945; Finney, 1945, 1946; Kempthorne, 1947, 1952; 
Cochran & Cox, 1950] and a considerable number of examples of frac- 
tionally replicated or confounded 2” or 3” factorial experiments are 
available for study. Few examples are to be found of factorial experi- 
ments in which the number of factor levels is a power of a prime. One 
example is given by Kempthorne & Tischer [1953] and is a 1/4 replicate 
of a8 X 4° X 2* split-split-plot factorial experiment. In the present 
paper a symmetric example of the 4” series is discussed. 


4° Experiment in Quarter Replicate 


The factors are denoted A, B, C, D and E, and the levels of each 
factor numbered 0, 1, 2 and 3. Corresponding with each factor we 
define a pair of pseudofactors each at two levels 0 and 1. Thus the 
four levels of each factor correspond with the pseudofactors thus: 


Level of A Levels of pseudofactors 
A’ A” 
0 0 0 
1 0 1 
2 1 0 
3 


and likewise with the other factors. 


Any treatment combination may be represented by the binary 
number, 
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where z; = 0 or 1. Successive pairs of digits correspond with different 
factors. Thus the treatment combination specified by 00, 01, 10, 11, 11 
represents the Oth level of factor A, Ist of B, 2nd of C and 4th of both 
D and E. 

The identity relationship chosen as the basis of a quarter replicate 
of this design is:— 


IT = A'B'C'D'E' = = (1) 


It is readily verified that no main effects or first order interactions 
are mutually confounded in this design. It must be remembered, 
however, that an interaction such as A’A”’B’B” although a third order 
interaction between pseudofactors is a first order interaction between 
factors. Successive multiplications of equation (1) yield:— 


A’ = BYC’D'E’ = = (2) 
A'A” = A"BIC'D'E’ = A'B"C"D"E" = (3) 
= A''C'D'E! = = BYC'C"D'D"E'E” (4) 
A’A"B'B" = A"B'C'D'E' = A'B'C"D"E" = C'C"D'D"E'E” (5) 
A'B! = C'D'E’ = A’A"B'B"D"E" = 


This shows that at most in (6) a first order interaction is confounded 
with a second order interaction. Similar equations hold for other 
pairs of factors owing to the symmetry of the design. Working out 
the details of the alias relationship in asymmetrical cases is tedious 
unless card sorting techniques are available [see Kempthorne & Tischer 
1953]. 

After defining a suitable alias relationship the treatment combina- 
tions admissible to the quarter replicate must be listed. In the present 
example it is not very time-consuming to carry out this by hand. 
Treatment combinations corresponding to each quarter replicate are 
specified by: 


Sum z; Sum z; 
Quarter replicate t= 1,3,5,7,9. | 7 = 2, 4, 6, 8, 10. 


First Even Even 
Second Even Odd 
- Third Odd Even 


Fourth Odd Odd 
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If it is decided to use the first quarter replicate defined above then 
treatment combination, 


01, 01, 01, 01, 11 
would be excluded since Sz; is odd and Sz; is odd. On the other hand 
combination 


00, 00, 00, 11, 11 


would be included since Sz; is even and Sz; is even. 

The first quarter replicate was used in a study of the effect of oestro- 
gen on the mitotic rate of the vaginal epithelium of the ovariectomized 
mouse (see Claringbold, 1956, for a full account of the work). The 
analysis of variance presented no unusual problems, having the form:— 


Source of variation Dé: 
Five main effects 5X3 = 15 
Ten first order interactions 10 X 9 = 90 
Remaining interactions as error 150 
Total 255 
DISCUSSION 


It is easy to see that the 4° factorial experiment is the smallest of 
the series which may be quarter replicated and yet leave interactions 
of less than three factors unconfounded: for confirmation check with 
equation (6) and ignore factor E, or alternatively consider the degrees 
of freedom (d.f.) in the 4* experiment. Main effects require 12 d.f. 
and first order interactions 54, a total of 66 which exceeds the number 
of observations in the quarter replicate, i.e. 64. 

The smallest experiment which may be 1/16th replicated is the 
4” with 1,024 treatment combinations. The 4° experiment may be 
quarter replicated using a similar alias relationship to (1) with the 
addition of F’, F’’ and F’F” to the second, third and fourth terms 
respectively. The 1,024 treatment combinations of both experiments, 
however, are probably too extensive for most experimenters. An 
admissible 1/16th replicate of the 4° experiment is not possible. It 
thus appears that the 4° experiment is the only one of the series which 
may be fractionally replicated unless further assumptions can be made 
about the main effects and first order interactions. 
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The three d.f. for each main effect may be reduced in many ways 
to sets of orthogonal single d.f. contrasts. One common method is 
to separate the sum of squares corresponding with the three d.f. into 
linear, quadratic and cubic components of regression. It may be 
possible to assume that the cubic component together with its inter- 
actions is negligible. For any factor at four levels we may compute 
the sum of all responses at each level. If these treatment totals are 


denoted T, to T, we may make the following orthogonal contrasts 
between them, 


The orthogonal contrasts giving the linear, quadratic and cubic com- 
ponents of regressions are respectively, 


L = 2A’+ A” = —37, —T,+T7; + 37, 

Q =A" 

C=2A” — A’ = + 37, — 37,+7,. 
It will be seen that the first set of contrasts is the set we would obtain 
if the four totals are considered the result of a 2’ experiment with the 
contrasts A and A” isolated as main effects and A’” as the interaction. 

If it is possible to regard the cubic component of regression as 

negligible, and if the contrast A’”’ is confounded then A’ gives an estimate 
of the linear component of regression and A’” of the quadratic com- 
ponent. It is thus possible to confound one of the main effect d.f. of 
‘the factor and its interactions with other factors. A design of this type 
is of use in optimum condition studies where the peak of a response 
surface is being sought, provided that the surface is no more’ complex 
than the second degree. With more severe assumptions of this type 


many more possibilities of fractional factorial experiments are prac- 
ticable iv. the 4” series. 
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SOME SMALL SAMPLE TESTS OF SIGNIFICANCE FOR A 
POISSON DISTRIBUTION 


C. RADHAKRISHNA Rao anv I. M. CHAKRAVARTI 
Indian Statistical Institute, Calcutta 


0. INTRODUCTION 


The large sample tests associated with Poisson distributions such as 
(i) goodness of fit, (ii) homogeneity of the observations, (i.e. arising 
from the same Poisson distribution), (iii) deviation in the frequency 
of zero etc., are not valid if the Poisson parameter is small unless the 
sample size is extremely large. In fact the chi-square approximation 
for the ‘index of dispersion’ (variance test for homogeneity) seems to 
depend largely on the magnitude of the Poisson parameter rather than 
on the sample size. Sukhatme’s (1938) model sampling results indicate 
that the sample size need not be large provided the Poisson parameter 
is not small. The values of the parameter chosen by him are, however, 
all greater than unity and it is not known what happens for smaller 
values. 

Fisher (1950) provided an exact treatment of the goodness of fit 
and variance tests by considering the conditional probability of the 
observations given their total, which in view of the sufficiency of the 
total is independent of the unknown Poisson parameter. The present 
paper is in a sense a follow-up of the technique given by Fisher. 


A number of problems have been investigated: 


(i) A test based on the likelihood ratio is suggested as an alternative 
to the variance test of judging the homogeneity of different observations 
or more specifically whether the observations arise from a single Poisson 
population and not a compound Poisson distribution. This test seems 
to be better suited for this purpose than the variance test in small 
samples*. 

(ii) An exact test of deviation in the ‘zero’ frequency or in general 
for any other frequency is given. The large sample test given by Cochran 
(1954) may not be valid for small values of the Poisson parameter when 


*It might be of interest to point out that the suggestion of using a likelihood instead of x? is quite 
old. [Fisher (1922), Neyman and Pearson (1928), Cochran (1936).] But the present application seems 
to be new. 
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the sample size is not very large. The exact test it may be noted is 
independent of the unknown parameter. 

(iii) Exact tests for goodness of fit and homogeneity have been 
worked out for the Truncated Poisson distribution. The appropriate 
index of dispersion test for the Truncated case has been given. This 
differs slightly from the treatment of David and Johnson (1952). 

Limited tables have been provided for the application of a number 
of tests associated with the Poisson distribution. 


1. HOMOGENEITY TESTS 
1.1. The index of dispersion 


If , are f observations from a Poisson population 
then it is known that given the total 


which is sufficient for the Poisson parameter, the conditional probability 
of the observations is (see Rao, 1952, pp. 36-37) 


al (1.1) 
which is multinomial with equal probabilities in f cells. To test the 
hypothesis that z, , x, , --- , 2, arise from Poisson populations with the 


same parameter we need only examine whether a total number of T 
observations from a multinomial distribution with equal probabilities 


could produce the frequencies z, , 2, , --: , Z;. The classical x’ test 
based on (f — 1) degrees of freedom for this is 
— E(z; | T)}? 
E(z, | T) 


which reduces to 


D= (x: — 


where ¢ = E(z,; | T) = T/f. This is referred to as the variance test. 
For the validity of the x’ approximation it is generally required that 
the expected frequency 7'/f should be greater than 5 but it appears 
that in the above case T/f could be much smaller than 5 as found by 
Sukhatme (1938) in his model sampling investigation. To investigate 
the nature of the approximation it is proposed to obtain the first four 
moments of the variance statistic subject to the condition that T and 
f are fixed. 
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The conditional expectation of any function ¢(z, , 22, --* , Zs) is 
connected with its total expectation by the following relation, where 
np is the common mean of the individual Poissons 


Hence 


E@|T) = {Coeff. of in eE(@)}. 


Therefore knowing the total expectation, the conditional expectation 
can be easily found. We can consider the statistic S = > xi — 
whose moments are known functions of » and then derive its conditional 
moments. For instance 


E(S|T) = {Coefficient ine*(f — 1)y} 


=2(8-4|7) = 4-0 


In the following equations expectation stands for conditional expecta- 
tion given 7’, which is dropped for convenience of notation. 
Similarly all the moments are evaluated: 


Function Exact expression x’-approximation 
wi(D) = 
= ag — (1 4) ay — 1) 
1 2 
8f 


+4y ny - 2 2) 
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[[48(f — 1) + 12(f — 1)’] 

(1 3) 48(f — 1) 

+ 48(f — 1) (1 2) 1)? 


+9(, _ 2) 


u(D) = 


8 (T+ (f — 6/2)’ 8 

__12 1 — (2f — 9)T + (f°/6 — 5f + 15) 12 

B(D) 8 = T(T — 1) f-1 

28, — 36, -6 = =3) — + 48f — 144) 


— 1) 


It may be noted that the expressions obtained above refer to x’-goodness 
of fit of a multinomial with equal probabilities in f cells and a total of 
T observations. For this case the expressions up to u, have been derived 
earlier by Haldane (1937) using a different method. Our object in 
reproducing these expressions is two-fold. Firstly, as an illustration 
of a simple method proposed for their derivation. Secondly, to provide 
a direct comparison of the exact moments and the #-functions with 
those of the x’-approximation. 

The rapidity with which the above 6-criterion tends to zero depends 
on the value of T'/f and to some extent on the size of f as well. Numer- 
ical evaluations of the criterion and the comparison of actual and 
approximate probabilities suggest that it may be misleading to use 
the x’-approximations if 7'/f is less than unity and that the approxi- 
mation yields good results for T/f over 3. Where exact probabilities 
are needed one has to use the multinomial probability (1.1) and cumulate 
the probability for values of D equal to or larger than the observed. 
Tables given at the end of the paper contain these expressions. 


1.2. Likelihood ratio test as an alternative 


If the observations arise from different Poisson populations then 
the likelihood of any set of parameters pu, , uo, --- , My IS 


II z,! 
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which attains the maximum 
-r 
IT 
while the maximum likelihood for a common uy is 


e(T/f)" + T] 
Hence the log likelihood ratio test for detecting heterogeneity is 


Ly = — > a, log, x; — T log, f + T log. T 


If T/f is large, —2L, is approximately x’ with (f — 1) degrees of 
freedom. According to general theory the statistics L, and D are 
asymptotically equivalent and there is not much to choose between 
them when 7'/f is large. Probably the same remarks about the size 
of T/f and validity of the x’ approximation as in the case of D hold 
for Ly also. But in small samples and/or when 7'/f is small indicating 
that an exact test is preferable, then the L, test seems to have certain 
advantages. The D statistic tends to be heavily grouped in small 
samples so that the D with a cumulative probability less than or equal 
to 5% may actually correspond to a much lower level of significance 
because of the gaps in D. With Ly, much closer percentages are ob- 
tained and consequently it has better chance of rejecting the null 
hypotheses. Ly, is preferable to D for discriminating a compound 
Poisson distribution from a simple Poisson. 

As we are using conditional tests the statistic may be simply defined 
by Li = > z; log, x, , just as in the case of D we may have D’ = >> 2’. 
Fisher (1950) gives the following observed distribution: 


Variate 0 1 2 3 Total 
Frequency 124 12 2 2 140 
The values of the statistics D’ and L are 
D' = Dorf; = 38 


Li = DX fii log, i) = 9.364362 


The cumulative probability using the exact expressions for D’ is .001112 
while for Ly it is .0005732 indicating that it is more sensitive. 


2. TEST FOR DEVIATION IN THE ‘ZERO’ FREQUENCY 


2.1. Large sample tests 


Situations arise in experimental investigations where the observa- 
tion ‘zero’ from a Poisson population is under-represented or over- 
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represented. When its frequency is unascertainable or when it is not 
considered we have the case of a truncated Poisson. It is therefore 
reasonable to examine whether the zero frequency is the offender when 
a goodness of fit test, such as x’, shows a departure of the observed 
distribution from what is expected on the basis of an underlying Poisson 
distribution. A natural approach to this problem is to decompose the 
total departure into two components, one due to the zero frequency 
and the other due to the rest of the frequencies. Both are important 
because of our interest in also examining whether the rest of the dis- 
tribution can be accounted for by a truncated Poisson. Departures may 
exist in both. This can be accomplished in large samples by a partition 
of the total x’*. If x} is the total x’? with x? for the truncated Poisson 
then the relation 


supplies the x? with one degree of freedom for detecting eignificams 
departures in the zero frequency. 

David and Johnson (1952) provide the following interesting dis- 
tribution of the number of decayed teeth. 


Number of decayed teeth O 12 3 4 5 6 7 8 91011 12 
Frequency of boys 61 47 43 35 28 15205 5 21 2 1 
Estimating the value of » for the Poisson by the formula 

m= T/f = 92 = 2.5736 


we find the value (after grouping the classes 7 to 12) 
= 152.3514 with 


As the classes 7 to 12 are grouped, the estimate of » used from the 
ungrouped data gives a slightly larger x’ for 6 degrees of freedom. 
As there is no proper method of adjusting for this an alternative estimate 
of » may be obtained by considering the classes over 6 as grouped [for a 
similar method in a contingency table, see Rao (1952) p. 199]. In the 
particular example considered above it is not likely to make any differ- 
ence because the classes grouped belong to the tail of a distribution with 
a small total probability. Further, the computations carried out by 
*In the case of a Poisson distribution the number of classes is theoretically infinite. But for pur- 
poses of computing x? it is necessary to group the classes beyond a certain value of the Poisson variable 
The large sample distribution of x* is valid when a finite number of classes obtained as stated above ir 
considered. In practice this is what is done. Since the tail frequencies will be small, they are combinec 


into a single class. In such a case the parameters have to be estimated from the grouped distribution 
if the x? approximation is to be valid. 
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Chernoff and Lehmann (1954) for a Poisson distribution truncated 
from above, show that no serious error is made. -Finally it -may be 
noted that x’ itself is an approximation and a small refinement in adjust- 
ment for grouping of classes with a-small probability is not worth 
attempting. 
Omitting f, = 61, u is estimated from the formula 
f — fo 


= 3.3431 


= 3.2075 
The goodness of fit x’ for the truncated case is 


x? = 31.0660 with 5d. 


Here again the same remarks as in the case of x3 hold good. The test 
for the ‘zero’ frequency is 


x = x2 — x? = 152.3514 — 31.0660 = 121.2854 . 


with 1‘d.f., which is extremely large indicating something obviously 
wrong. But it must also be noted that even the hypothesis of a trun- 
cated Poisson is not tenable, as x’ is also significant. Probably this 
is a case of a compound Poisson with further trouble at the ‘zero’. 
The hypothesis of homogeneity of the observations from a truncated 
Poisson is examined in Section 4. 

Alternatively we could use the likelihood ratio goodness of fit test 
and iis decomposition. If fy , fi , --- denote the observed frequencies 
and ~- ™, - - those estimated under a given hypothesis then 


Le(p) = log. 
r>0 r 
For the sruncated Poisson 


Lal) = log, 


‘where m;’s are estimated frequencies for the truncated Poisson. The 
difference L, = Le(p) — Le(t) provides the test for zero frequency. 
The statistic —2L, is asymptotically x’ with one degree of freedom. — 

One could use the large sample test for an individual frequency 
suggested by Cochran (1954) or the equivalent likelihood ratio criterion 
for examining the zero frequency directly but such a test will not reveal 
the true situation unless the rest of the distribution is known to behave 
as a Poisson distribution. It is with this practical purpose in view that 
the statistics x? and L, are defined. 
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2.2. Exact treatment of the ‘zero’ frequency 


In Section 1 it was seen (Formula 1.1) that the conditional probability 
of the Poisson variables is multinomial with number of cells equal to 
the number of observations (f) and sample size equal to the sum of 
the observations (7). The zero observation corresponds to an unre- 
presented cell in a sample of size T. The problem can be posed in a 
familiar way as one of the classical occupancy problems mentioned in 
several text books on probability. From an urn containing equal 
numbers of balls of each of f colours, T balls are drawn. What is the 
probability that f — r colours are unrepresented? This probability is 
easily found to be 


r=0,1,2,--- 


Or using the notation of differences of zero the above probability can be 


written 
P, = (i Jaro" 
rw 
If T < f, Pr.; automatically takes the value zero fori = 1, 2, --- . 


To test the hypothesis that the observed frequency f, is more than 
the expected we calculate the probability / 


I-fe f 
> 
1 I-fe-1 
using the convenient expression. If this is less than the level of signifi- 
cance chosen we reject the hypothesis. To test the hypothesis of under- 
representation of zero the probability to be calculated is 


I-fo-1 f 
1 
When both positive and negative departures of f, from the expected 
are considered there is some arbitrariness in judging significance. 
One may follow the rule of rejecting any observed f, not falling in the 
range (0, , 0.) where 0, and 0, are determined from the formulae 


f-0,-1 a a 
1 f-0,+1 


in which a stands for the level of significance. 
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The values P, can be obtained from the comprehensive tables for 
A'0” = r! prepared by Stevens (1937) and reproduced in Fisher and 
Yates’ Tables. Stevens considers several examples from genetics where 
this test is useful. In Fisher’s example we find using Stevens’ tables 
for T = 22, 


16 


> P, = .00158237 
which shows that the ‘zero’ is somewhat over-represented. It is not 
easy to compute the expressions for cumulative probabilities for large 
values of 7. Normal approximation can be used if f is large*. The 
exact expressions for mean and variance are found below. 
The k™ uncorrected moment of r is given by the value of the follow- 
ing expression at z = 1 


ss — (2f - 4)’ 


K(f.) - (5) 
For large f, the normal deviate is 
{re - + 


In David and Johnson’s (1952) example, the frequency for zero is 
fe = 61 out of a total of f = 265. The value of T = >> f,z,; = 682, 
= 2.5736 


61 — 201.095 


i= = 


VV Go) V146.800355 


*One of the referees of this paper points out that the fact that normal approximation can be used 
for f large is proved by I. Weiss in 1952 in a thesis entitled ‘‘Limiting distributions in some occupancy 
problems"’. (Technical Report No. 28 prepared under tract N6onr-25140 for the Office of Naval 
Research at Stanford University, July 1955). 
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which is too large for a normal deviate, establishing significance in the 
frequency of zero. 


If T is large compared to f then P, can be approximated by [see 


Feller, (1950)]} 


where \ = f e~””’. In this case the Poisson tables may be used for 
obtaining the cumulative probabilities. 


3. CONDITIONAL VARIANCES AND COVARIANCES OF THE FREQUENCIES 


Any other frequency or in general any linear combination of the 
frequencies can be tested for deviation from the expected values by 
using the exact variances and covariances and the large sample normal 
approximation. By using the technique given in Section 1 the condi- 
tional mean and variance of f; , the frequency of the observation 7, is 


From the above formulae the exact mean and variance for any linear 
function of the frequencies can be calculated and its significance can 
be tested by using a normal approximation if f is large. In an example 
given by Cochran (1954) where f = 240, T = 388, the value of f; = 52. 
To test whether it significantly departs from the expected we compute 


= _ 52 = 33.60 _ 
VVG;)  ~~23.8144 


which is significant. The large sample approximation given by Cochran 
(1954) slightly over-emphasizes significance. The exact probability 
for f; quoted by Feller (1950) is 


the summation extending over those j > f; for which j < f and ij < T. 
This is, however, difficult to deal with computationally. 
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4. TRUNCATED POISSON 


4.1. The goodness of fit test 


In the case of the truncated Poisson the probability of the observed 
frequencies f, , fz, --- (with a total of f’) given T = >>; f; is 
T! 1 
f) 


where 


= 4/0" 


o(T, f’) 


The corresponding expression for the full Poisson is simply f”. The 
probability for f,; , f2 , --- given T is independent of the unknown 
parameter because of the sufficiency property of the statistic 7 even 
for the truncated case. The likelihood ratio statistic or the x’ statistic 


where 


l-—e”r! 


and m is estimated from the formula 


m, = 


f 


can be used to test goodness of fit in large samples. The remarks made 
in Section 2.1 about the validity of x’ approximation hold good in 
this case also. In the case of the likelihood ratio statistic it has to be 
multiplied by (—2) before referring it to the x’-table. If the cell 
frequencies are small the exact cumulated probability for any one of 
the statistics can be calculated by using the formula for exact prob- 
abilities of the different configurations. With small samples the likeli- 
hood ratio test is preferable. 


4.2. Homogeneity tests 


The likelihood ratio test for homogeneity based on the observations 
Zi , from f’ truncated Poisson populations is 


a 
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Ly(t) = log. mi 


where 


mM; 


The statistic —2L,,(¢) can be used as x’ with (f’ — 1) degrees of freedom 
in large samples. In small samples the exact probabilities can be 
calculated from the conditional probability 


T! 1 
where ¢(7', f’) is as defined in section (4.1). 
The analogue of the index of dispersion statistic in the truncated 
case appears to be 


(2; — #)? + + m— 2) 


where = m/(1 —e”). In large samples this has a x’ distribution* 
with (f’ — 1) degrees of freedom. In small samples the statistic to be 


used can be simply > 2° . 


Using the illustration of David and Johnson (1952) we found the 
estimate of u to be 3.2075. The index of dispersion D(#) for the trun- 
cated case is 


— + + m— = 1009-9804 _ 349 5952 


based on 203 d.f. Using the normal approximation for x’, the normal 
deviate is found to be 5.27. This is somewhat higher than the value 
4.95 of the statistic suggested by David & Johnson, which differs from 
D(® in the denominator being simply Z. 
4.3. Exact moments of x 
In the truncated case, mean and variance of >> zi are 
ED x) = — — 2, f) + o(T — 2, f’ — 
+ 
wa = (T (f(T — 4, f’) + — — 4, f’ — 1) 


*This follows from the general theory of the x? test discussed by Rao (1952, p. 177). 
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+ + — 3, f’) + 2f'(2f’ + No(T — 3, f’ — 1) 
+ 2f'(2f’ — Ie(T — 3, f’ — 2)} 
+ + — 2, f’) + f(2f’ + — 2, f’ — 1) 
9) 
+ — 1, f) —1,f' — 
— 

where 


= —1)---(T—r+})) 


5. SOME EXACT TESTS FOR THE BINOMIAL DISTRIBUTION 


It may be noted that investigations similar to these in the Poisson 
case can be carried out for the Binomial, Normal and other distributions 
where the nuisance parameters can be avoided by considering conditional 
probabilities. If f sets of s trials are made with probability p of success 
in a trial then the probability of frequencies f, , f: , --- f, of the possible 


Observing that >> rf, = T is sufficient for p and has the probability 


Ss! 
TS — 


we find the conditional probability of the frequencies 
f! — T)! (*) 
fol --- S! I] 


Using this expression the exact probability of any criterion of goodness 
of fit can be computed. 
The probability of f. given T and f is found to be 


S= sf 


This expression can be used in computing the probability of tests of 
significance concerning the frequency of zero success. 

If f is not small, normal approximation holds good for fy . The 
exact mean and variance of fo are 
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(S-7! 


— — — 2s)! 
V(fe) = E(fo) — (EGP + — Ti 


6. TABLES AND THEIR USES 
6.1. Critical values 


Critical values for a level of significance a < .05 of the statistics 
> 2? , } z, log, z, , >> f, log. (f,r!) and frequency of the ‘zero’ class 
are given for T = 3(1) 10 and f = 3(1) 10(10) 100. The exact level of 
significance occurs as a lower entry under each critical value. When 
the probability levels are too low, the next lower value of the statistic 
is also recorded with the corresponding probability level if it is not 
much above 5 percent. If the observed value of a statistic happens to 
be equal to or greater than the tabulated value, the hypothesis is to be 
rejected at a level of significance a < .05. 

It will be seen that in many cases the exact probability (level of 
significance) is much below 5 percent. This is due to the fact that the 
distributions of the statistics under consideration are heavily grouped; 
the inclusion of a smaller realisable value of the statistic would increase 
the probability beyond the 5 percent point. 

The statistic >> f, log, (f,r!) is comparable to x” when the expression 
(T — f) log. f + T — T log, T is added and multiplied by 2. 

For intermediate values of f, critical values of the statistics may be 
obtained approximately by choosing the nearest f tabulated values 
or by interpolation between two values of f. 


6.2. Illustrative examples 


The following frequency distribution is used to illustrate successively 
the use of Tables 4, 1, 2 and 3. 


r 0 1 2 3 
I. 63 5 1 1 
Here f = 70 and T = > rf, = 10 


(i) Test for ‘zero’ frequency (one sided test to examine an excess in 
the expected zero frequency): Observed frequency of the zero 
class fp. = 63. The tabulated value for 7 = 10 and f = 70 is seen 
to be 63 with the associated significance level of .01. Significance 
is noted when the observed value is equal to or greater than the 
tabulated value. In this case the deviation from the expected is 
significant indicating an excess in the expected frequency of zero. 
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(ii) Tests for homogeneity: Two statistics }> 2; and >> z, log, z, are 
computed to test whether the individual observations come from 
the same Poisson distribution. 


(a) = =154+414+9.1 


= 18 
This is in excess of the corresponding tabulated value 16 for 


T = 10 and f = 70 with a significance level of .03. So the hy- 
pothesis of homogeneity is rejected. 


(b) > 2; log. z; = Dorf, log.r 
4.345 


The tabulated value of >> 2; log, z; for 7 = 10 and f = 70 is 
3.2 at a significance level of .03. So the hypothesis of homogeneity 
is rejected on the basis of the second test. 
(iii) Test for goodness of fit of the Poisson distribution: 

The statistic >> f, log. (f-r!) is computed to test goodness of 
fit and this comes out to be 271.549. On referring to Table 3 
for the critical values of this statistic, it is seen that the observed 
value is in excess of the tabulated value of 271.29 at a level of 
significance equal to .03. So this indicates that the observations 
do not belong to a Poisson distribution. 
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MATCHED PAIRS IN SEQUENTIAL TRIALS FOR SIGNIFI- 
CANCE OF A DIFFERENCE BETWEEN PROPORTIONS 


W. Z. Bittewicz 
Obstetric Medicine Research Unit (Medical Research Council), 
Midwifery Department, University of Aberdeen, Scotland 


INTRODUCTION 


I. Bross (1952) and P. Armitage (1954) have discussed the applica- 
tion of sequential methods to medical trials, pointing out that not only 
are they economical in use but also keep the effectiveness of a given 
method of treatment under continuous review. Yet sequential methods, 
though firmly 6stablished in various other experimental fields, have 
not often been used in medical research. 

One of the reasons for the limited use of these methods is that the 
success of the treatment and its length are frequently related to such 
individual characteristics of the patient as age, sex and severity of 
condition at entering the trial. The chance of contracting a particular 
disease may also vary with such environmental conditions as housing, 
dietary habits and occupation. For such reasons, medical research 
workers may insist on having experimental and control cases similar 
with respect to characteristics which they regard as important in a 
given case. From the statistician’s point of view this approach implies 
single or multiple stratification and a certain degree of correlation in 
response to treatment between members of the same pair. 

Stratification means a more complicated sampling machinery than 
would be necessary for unrestricted pairing, and correlation raises the 
question of the applicability of the usual sequential formulae to a 
matched-pair trial. The usual test based on unrestricted pairing 
will show, with required probabilities, whether the new treatment is 
or is not significantly better (or worse) than the standard, without 
involving us in any difficulties of selection or theoretical problems. 
At first glance, therefore, it might seem that the introduction of matched 
pairs is an unnecessary complication unless it can be shown to result 
in a sizeable reduction in the average amount of testing required to 
reach a decision. 

Even if there is little or no reduction in the average sample size, 
the matching procedure may be worth using for the following reasons. 
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When a research worker insists on matching he usually has two things in 
mind. The first is to eliminate as far as possible those differences 
between the subjects which may confuse the issue, so that when the 
experiment is completed the difference between experimental and control 
groups can be attributed to differential effect of treatments alone. The 
second, less frequently stated, is the possibility that the response to 
the new treatment may vary between strata and that some information 
on this point may help to assess the results or bring to light some 
hitherto unsuspected characteristics of the new treatment. 

On occasions the statistician may suspect that the research worker’s 
request for matching really stems from fear that some ‘‘unusual pairing”’ 
will vitiate the results and will consequently reassure him on this point. 
When, however, there are clear cut differences between strata it seems 
wasteful not to use this additional information. In addition the clarity 
of the experiment is obviously improved when extraneous factors 
known to “affect’’ the issue are controlled. 

The usual experiment based on unrestricted pairing can shed no 

light on the differential response of strata while a careful scrutiny 
of a fairly long matched-pair experiment will yield some information 
on this point. To sum up, when there are good grounds for it, matching 
results in an improvement of the experimental design both from the 
medical and from the statistical point of view; it increases the confidence 
with which the research worker will interpret the results; it uses all 
information at our disposal; and the analysis of a completed experiment 
may yield valuable information, otherwise unobtainable. Thus, the 
necessity of setting up a slightly more complicated sampling machinery 
does not seem to be an excessive price to pay for these advantages. It 
will be shown below that there is no additional cost in terms of sample 
size. 
The usual formulae for the sequential test for’ significance of a 
difference between proportions are valid when there is no correlation 
in response between members of the same pair. The object of this note 
is to examine the case when, as in matching, such a correlation exists. 
The note arose from a practical situation in which a sequential trial 
was strongly indicated on economic grounds, while the similarity of 
experimental and control cases was regarded as very important from 
the medical point of view. 


DEFINITIONS 


For convenience, open, single-tail plans will be considered here. 
The procedure can be easily extended to cover two-tail tests. 
Given the critical proportions of successes, p and p*, (p < p’*) for 
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the standard and the experimental treatment and the acceptable errors 
of both kinds*, a and 8, the sequential plan is fully determined. The 
probabilities of success determine the value of 


which together with errors of both kinds determine the upper boundary 
line. For the case of ‘no difference’, p = p* is accepted here so that 
for the lower boundary 6, = }. 

In a sequential trial subjects are selected in pairs, one member 
being allocated to the experimental and the other to the control group. 
Of all pairs only the untied pairs yield information about the difference 
between the experimental and the control groups. Thus only SF and 
FS pairs contribute to the test, where S stands for “success”, F for 
“failure” and the first letter refers to the control group. It will be noted 
that @, is the expected proportion of pairs favourable to the new treatment 
(FS) among all untied pairs if the population values of the proportions 
of successes in the two groups x and x” are equal to p and p’ respectively. 

Using the graphical representation suggested by Armitage (see 
Figure 1) the equationg of the boundaries are given by: 


Upper boundary y, = a, + sz 


(2) 
Lower boundary y, = —a, + sz 
where z is the number of untied pairs. 
When 3) 
1 0; 
a, = (4) 
logio 1 6; 
OB10 46,(1 6;) (5) 


*The adopted values of a and § depend on the appreciation of risk involved in coming to a ‘‘wrong” 
decision. a is the maximum proportion of instances in a long series of tests in which we may be led to 
conclude that the experimental treatment is superior when this is not true. § is the maximum proportion 
of instances in a long series of tests in which we may be led. to conclude that the experimental treatment 
is not superior when in fact the opposite is true. In other words pairs of treatments for which the true 
@ < 4 will lead on the average to a decision in favour of the standard treatment with relative frequency 
of at least 1 — a, while pairs of treatments far which the true @ > @: will lead on the average to a 
decision in favour of the experimental treatment with a relative frequency of at least 1 — 6. 
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where 
l-—a 


B 


t = logic z = logio 


a and £6 being the adopted errors of both kinds. 

In the usual method pairs are drawn at random from the whole 
population. When matching is used pairs are still drawn at random but 
both members of each pair belong to the same randomly selected 
stratum of the population. 

If the population is divided into R strata and pairs are drawn from 
them at random the probability of drawing an FS pair is given by 

* a:pi(1 — p;) where p? and 7p; are the critical proportions of suc- 
cesses in the strata and a,’s their proportional sizes. The probability 
of drawing an SF pair is given by >> a:p;(1 — p?). Consequently in 
an experiment with matched pairs the critical value of 64, is 


a,pi(1 
a(pi(1 — + pl — pi) 


Ou, = 


where 
R R R 
aD; =?D, aD; = and a; = 1 (6) 


When 7; = 7, in all strata, that is when there is no difference between 
the experimental and the control groups, 64, = } and for the overall 
proportions of successes we have p = p’. 

Given the value of 6, the boundaries for a matched-pair sequential 
plan are calculated from expressions (2)—(5). 

The average number of untied pairs vi, needed to reach a decision 
depends on the value of @ in the population. Three points are of special 
interest and may be estimated from the expressions below. For an 
experiment with unrestricted pairing we have: 

When the true value of @ = @, 


a, — B(a, + a) 
20, —s-1 7) 


When the true value of 6 = @, = 4 
a, — a(a; + 


Ne, = 


= (8) 
At maximum, when @ about half-way between } and @, 
= (9) 


R 
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When the true value of @ is greater than 6, or smaller than } the average 
fi, is smaller than fis, or vip, respectively. For a matched-pair experiment 
the corresponding values of , and are obtained by 
substituting in the above formulae the corresponding values of param- 
eters of the matched-pair plan. 

To estimate the total amount of testing, that is the total number 
of tied and untied pairs required to reach a decision, we multiply the 
appropriate value 7, by Ky . 

We have thus for the trial with unrestricted pairing: 


Ms, fe, Ko, = he. Ko, and Nous = Tomax x Kons: (10) 


For a matched-pair experiment substitute the corresponding values for 
Ox 


In the case of the usual experiment 


1 
(11) 


For a matched-pair experiment we have 


K, = 


Ku =- (12) 


The probabilities used to estimate K are not necessarily the critical 
values of p and p* but the values of x and x* in the population having 
@ given value of true 6 or 64% . Of course when the population @ is 
assumed to be equal to 6, as in (7) the values of x and x” are the same 
as the critical values used to determine the boundaries. 


SEQUENTIAL PLANS FOR MATCHED PAIRS 


When there are sound reasons for having pairs similar in certain 
characteristics we may assume that there is enough information to 
enable us to obtain estimates of a; and p; reasonably close to the popu- 
lation values. The proportion of successes thought to indicate superi- 
ority of the new treatment will usually be given in terms of p’*, while 
the method requires knowledge of the values of p; for each stratum 
separately. Two of the possible methods of determining these values 
will be considered: 

(1) We may assume a constant value of @ in all strata, i.e. a par- 

ticular set of pj values. 

(2) We may assume a proportional change in the proportion of 

successes or failures in all strata, i.e. a particular set of @; values. 
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Method 1. 


The assumption of constant @ in all strata implies that the greater 
the proportion of successes already obtained in a stratum, the smaller 
the proportional improvement necessary to make the new treatment 
worth while. The absolute difference between p; and p,; depends on 
the value of p; and is generally greatest when p, is not far from 4. 

It will be shown later* that for any set of strata and any set of p; 
values 64, > 6, where @y, is defined by (6), and also that for any set 
of true proportions of successes, x; and x; , resulting in a given value 
of true 6, constant in all strata, the average number of untied pairs 
and the average number of all pairs required to reach a decision is in 
general smaller for the matched-pair trial than for the usual method. 

The value of constant 6,4, may be found from the equation of the 
type: 


The solution of this equation is rather cumbersome and it is much easier 
to proceed by successive approximation as shown in Table 1 below. 


TABLE 1 
Calculation of Constant @4,, 
Strata 
Item All 
I II Ill strata 
Proportional sizes of the strata, a; 2 6 3 1.0 
Standard treatment, p .20 .50 .70 .50 
New treatment critical value, p* — — —_— .70 
Impliedt values of p7 for 
Ou, = .72 .3913 .7200 .8571 . 6954 
Ou, = .73 .4088 . 7300 8632 7057 
Ou, = .725 .3973 . 7250 8602 .7000 
{These values are obtained from 
Ou Pi 
1 — p; — Ou,(1 — 2p,) 


*Bee section: Proofs. 
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Thus for the critical value of p* = .7-this particular set of strata sizes 
and p, values leads to 6y, =. .725. Arriving at the value of Oy, by 
successive approximation rather than by solving equation (13) has 
the advantage of showing at the same time the values of p; in the 
strata. Assuming the acceptable errors to be a = .025 and B = .05 
and substituting the value of 6, = .725 into (3)—(5) we obtain: a, = 7.5, 
a, = 6.1, and s = .2334. The boundaries of the plan are then given . 
by (2), as: 

7.5 + .23342 
and 


Y, = —6.1 + .23342 


The sequential chart corresponding to this plan is shown in Figure 1. 


SEQUENTIAL CHART 
+2334 
184 42-025 
vw 


30 35 40 
NUMBER OF UNTIED PAIRS, x 


NUMBER OF FS PAIRS MINUS NUMBER OF SF PAIRS, Y¥ 


FIGURE 1. SEQUENTIAL CHART FOR CONSTANT @ = .725; SHOWING ONE OF THE 
EXPERIMENTAL SAMPLE PATHS. 


Only the untied pairs are plotted on the chart. The sample path 
starts at the origin and for each pair favourable to the new treatment 
(F'S) ascends by one step, while for each unfavourable pair (SF) it moves 
by one step downwards. Thus the first three untied pairs of the sample 
shown in Figure 1 were FS pairs, the fourth was an SF pair, and so on. 
Sampling proceeds until the sample path crosses one of the two bound- 
aries. If the sample path crosses Y, we may conclude that the true 
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difference between the treatments ig not zero; whereas if Y, is hit we 
may conclude that the difference, if any, is smaller than the critical 
value chosen at the beginning of the trial. 


Method 2. 


If we want to base our plan on proportional changes in the strata 
we may take either the proportion of successes or the proportion of 
failures as the starting point. If we assume that the proportion of 
successes changes proportionately in all strata we have 


Pi Pr Pp 

This implies that the greater the proportion of successes already observed 
in a stratum the greater is the absolute improvement necessary to make 
the new treatment worthwhile. It may, of course, happen that a 
given p*/p = will be found impossible for some strata or will lead to 
unreasonable values of p; . However, the problem may always be set 
in terms of a proportional decrease of failures, so that 


1— pi 


In this case, however, the implications are different. Proportional 
decrease of failures implies that the greater the number of successes 
in a given stratum the smaller the absolufe improvement required to 
accept the new treatment. In this respect the approach resembles that 
of Method 1, in which the proportional improvement decreased with 
the increase of the proportion of successes alregdy observed in a given 
stratum. To construct a sequential plan on the assumption of pro- 
portional changes we determine @y, , from (6) where pj; = Ap; and 
= p’/p. If proportional decrease of failures,is assumed, 04, is ob- 
tained by substituting in (6) values of 

The remaining steps are the same as described for Method 1. It can 
be proved that 8, > 6, and that the average sample length is smaller 
than for unrestricted pairing. 


SAMPLING EXPERIMENT 


Ten samples were drawn, with the aid of tables of random numbers, 
according to the rules of each method, to see how the two methods 
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worked under sampling conditions and how they compared with the 
usual method of unrestricted pairing from the whole population. 

The sampling population and changes in the proportions of successes 
implied by each method are shown in Table 2. 


TABLE 2 
Experimental Population and the Assumed Changes in the 
Proportions of Successes 
Strata of the population 
Item All 
I II Ill strata 
Sizes of the strata, a; 2 5 3 1.0 
Standard treatment, p .20 .50 -70 .50 
New treatment, p* 
Unrestricted pairing — .70 
Method 1, constant 
Ou, = .725 .367 .70 
Method 2, proportional decrease 
of failures \’ = 1.666 .52 .70 .82 -70 


If a fixed-size sample were drawn, then with errors a = .025 and 
8 = .05, 150 pairs would have to be drawn, yielding on the average 
75 untied pairs. 

Details of the sequential plans are set out in Table 3. 

The results of the experiments are shown in Table 4 where samples 
are arranged in order of length. The average length of the three groups 
of samples agrees very well with the expected values. There appears 
to be nothing unusual about the scatter of individual sample lengths. 
The average number of untied pairs was 43.3, 31.8 and 35.4 for un- 
restricted pairing, constant 6, and the proportional change methods 
respectively; all three values agree closely with expectations given in 
Table 3. A check on the average proportions of SS, FF, FS and SF 
pairs showed that these too were in good agreement with theory. 

In any individual sample the proportions of cases drawn from 
various strata may differ a great deal from those assumed. Examination 
of experimental samples did not give any clear evidence that the usual 
sampling variation in the proportions of cases drawn from variaus 
strata was associated with the length of the sample. In the case of 
samples drawn on the basis of proportional change there was, however,:a 
suggestion that samples in which Stratum III, with the lowest 6; and 
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TABLE 3 
Details of the Sequential Plans 
Unrestricted Proportional 
Item pairing Constant 04 change 

6, or Ou, .700 .725 .716 

a 8.6 7.5 79 

a 7.0 6.1 6.4 

8 . 2058 2334 . 2234 
Average number of untied pairs 

When 6 = 6, or 04 = Oy, 40 32 35 

When @ = 0y =} 32 25 27 

At maximum when 6 or 6a, 

approximately half way between 

4 and 6 or Oy, 63 48 53 
Multiplier K or K y* 

When 6 = 6, or 64 = Ow, 2.00 2.25 2.16 

When 6 = 64 = 4 2.00 2.27 2.27 

At maximum 2.00 2.27 2.21 
Average number of all pairs 

When 6 = 6 or 064 = Oy, 80 72 75 

When @ = 0y =} 64 56 61 

At maximum 126 110 118 


*For the purpose of calculating K, pi values were treated as the true proportions of successes as 
indeed they are arranged to be in this experiment. The equality of K values in the first column is 
accidental; when = .5 we have K = 2 irrespective of the value of r*. 


highest K, was over-represented tended to be slightly longer. One 
would expect Method 1 to be the more stable of the two, since though 
K varies from stratum to stratum the critical ratio of FS to all untied 
pairs remains constant. In Method 2 both K and 6 vary from stratum 
to stratum and in extreme strata the value of the latter may differ a 
great deal from that used to determine the boundary lines. 

To check this point five additional samples were drawn by each of 
the two methods; drawing from the strata with the same values of 
p; and p; , but with probabilities .2, .3 and .5 instead of those given in 
Table 2. The average length of the five samples in case of constant 
6, was 76.4 pairs (expected length 72 pairs); the shortest series was 
54 pairs while the longest was 108. In the case of proportional change, 
however, the five samples yielded an average of 144.4 pairs (expected 
length 75 pairs) and the shortest series was 92 while the longest took 
203 pairs before a decision could be reached. This result confirms the 
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inference that Method 2 is more sensitive than Method | to the variation 
in the proportions of cases drawn from various strata. The importance 
of variation in the proportions of cases drawn from various strata 
should not be exaggerated since the proportions were deliberately 
arranged to do most harm. Stratum III, i.e. the one with the highest 
K in case of constant 6, and with the highest K and the lowest 6; in 
case of proportional change, was heavily over-represented. This par- 
ticular distribution of cases served to make the point clear but in practice 
such extreme variation occurs rarely and moreover would be very 
unlikely to persist as the test goes on. It appears that under normal 
sampling conditions both methods of arranging a sequential plan for 
matched pairs yield consistent results, the method of constant 6, 
being slightly more stable. 


PROOFS 


Various statements made in the preceding sections will now be proved. 

It will be shown now that 64, > 6,. Oy, (6) was defined in such a 
way that >>" a,p: = p’, the critical value for the new treatment. The 
values of p and p* determine 6, for an experiment with unrestricted 
pairing. We will show that, if 6, were assumed to operate in all strata, 
the resulting >°1 a,p; would in general be smaller than the critical 
value p* for which the plan is to be designed. It follows from (1) that 
the difference between the critical value p* and >* a:p; = pix may be 
written in this case as: 


Pu 1 —p— 6(1 — 2p) 
Remembering that p = )>-* a,p, we arrive at 


6(1 — 6,)(20, — 1) 


p Pu = — 2p) 
— pi) 


the summation its taken over all 7,7 combinations. Since } < 6, < 1, 
this expression must be positive unless all p; = p when it becomes zero. 
It follows therefore that the constant @,, used for a matched-pair plan 
must in general be greater than @, and that the greater the differences 
between strata the greater the difference between @y, and @, . 

To prove the same proposition for Method 2, that is assuming 
proportional change in 7; in all strata, it is easier to start directly with 
the difference 04, — . 
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If the proportionality factor is \ > 1, we may write from (1) and (6) 


p> adp(1 — — p) 


R ~ 1+A— 
a;[Ap(1 — pi) + — Ap,)] 


1 
which leads to 


Ou, aad 6; = 


AA — 1) awa,(p; — 
(15) 
(1 + A — 2dp) > — pi) + — Ap,)] 


1 


which is positive unless all p; = p, and therefore 04, > 6, .* 

Substituting the true @ of the population for 0, , and zm for p, we 
may deduce the relation between the true @ and the true @y on the 
assumption that the proportion of successes in the strata changes 
proportionately to */r = X or according to the constant 6 rule. 
From (14) and (15) we see that if @ = 4 (ie. \ = 1) then whatever the 
values of 7; , 94% = 3; when 6 > 3 then 0y > 6, and when @ < 3} we 
have 64 < 6. These relations will be found helpful when we pass on 
to consider the average sample size. 

For any sequential trial the average number of untied pairs iis 
required to reach a decision depends only on the true value of 6. Three 
points of interest on the Average Sample Number curve (untied pairs) 
can be obtained directly from the parameters of the plan; the formulae 
are given as (7), (8) and (9). Using (3), (4) and (5) we may write these 
formulae for 7i, as decreasing functions of 6, or Oy, . 

Therefore for the three points at least the average number of untied 
pairs needed to reach a decision is smaller for a matched-pair plan than 
for an experiment with unrestricted pairing since 04, > 6, . Other 
points of the Average Sample Number curve, except the two limiting 
cases of 6 = 0 and 6 = 1, cannot be obtained directly from the param- 
eters of the plan but have to be calculated indirectly. 

Figure 2 shows the Average Sample Number curves (untied pairs) for 
the experimental plans. The scale of 6, has been so arranged as to 
show the values of 6,4 and of the average sample size corresponding 
to given values of @ for the experiment using unrestricted pairing. 
The Average Sample Number curve for the constant @y plan is pulled 
slightly out of shape as the result of this arrangement. The two curves 
show that for any value of @ in the population and its corresponding 


*The same holds when a proportional decrease of failures is assumed. All that is required for the 
proof is simultaneous switching around of the designations ‘‘success,’’ ‘‘failure’’ and ‘‘control,”’ ‘‘experi- 
mental.”” Such a change of ‘‘names” does not affect @ so that the problem can then be re-stated with the 
usual designations. 
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AVERAGE SAMPLE NUMBER CUNTIED PAIRS) 


UNRESTRICTED PAIRING 


CONSTANT @ 


NUMBER OF UNTIED PAIRS 


2 3 4 4 7 
*2 3 “4 7 9 


‘725 


FIGURE 2. COMPARISON OF ASN CURVES FOR SEQUENTIAL PLANS SPECIFIED IN 
TABLE 3. 


value of 64 , the experiment using matched pairs is likely to be the more 
economical. A similar picture can be obtained for the proportional 
change method. 

The average length of the whole experiment depends on the true 
value of @ in the population and on the values of x and x” for the two 
treatments. Since by matching we make the numbers of each pair 
similar in response the untied pairs must become less frequent; it is 
therefore not surprising that it can be shown that K,, > Ky unless 
all az; = 

For matched-pair plans, part of the gain in the average number of 
untied pairs is offset through K,, > Ky ; we have therefore to consider 
the total sample size. 

Suppose x; = p,; and x; = p; where x; and x; are linked by a constant 
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64, ; in other other words suppose that our estimates of proportions 
of success in the strata and our guess of the efficiency of the new treat- 
ment are correct. We want to know whether the matched-pair plan 
is more or less economical than the usual method based on unrestricted 
pairing. Let us enquire for what kind of population structure the average 
total sample size is at maximum for a matched-pair plan. 

From (7) and (12) by substituting (3)—(6) we may write* 


1 1| api(1 — pi) 
Q= api(l — pi) logio = 
ap.(1 — pi) 
- logio — 
2 — pi) 
where 


C= logio — B logio i= 


a and 6 being the errors of both kinds adopted for the trial. 


For the maximum of Ny, = Key jfie,, we want Q to be a minimum. 
Write 
R 
Z, = = pi) 
(17) 
Z, = api(l — 
then z, — z, = p* — p = 2u, a constant for any given problem. 
Consequently 
1 Z. Z Z. 
Q= C logio Z, — (Z, + logio 
1 (18) 
= C [Z; logio + (Z, + 2u) logio (Z, + 2y) 


— 2(Z, + logs (Z: + »)] 


*This proof of the propositien was suggested to me by Dr. M. R. Sampford, A. R. C. Unit of 
Statistics, University of Aberdeen. 
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Differentiating with respect to Z, , we find Q to be a decreasing function 


of Z, , so that Q is a minimum when Z, is a maximum. We may write 
Z, as 


Z, = p — pp — — (19) 


For maximum Z, we want the expression in the bracket of (19) to be 
a minimum. Since this expression is zero when all p; = p it remains 
to be shown that it is positive for p; ¥ p. 

It was shown by other means (14), that 04, > 96, ; let us write 
this in terms of Z’s, 


Z2 p(l — p) 
p(l—p)+p(l — 


~ + — p=) + — 


Since the denominator of this fraction is always positive the numerator 
is also positive for p; * p since 6y, > 6, . Consequently the bracket 
of (19) is always positive for p; ~ p. Therefore Z, is a maximum when 
all p; = p. Thus the average total sample size for a matched-pair 
plan is at maximum when p; = Pp and it is then equal to the average 
total sample size of a plan with unrestricted pairing. 

When the experimental treatment is just as good as the standard, 
Le. p; = 7; = writex = 1 — 46,(1 — and y = 1— 40y, (1 — Oy,). 
Then from (8), (11) and (12) 


= p(1 — p) log, (1 — x) — ap(1 — p) log,(l1—y) (21) 


Ou, — 


(20) 


By expanding the logarithmic series the left hand side can be shown to 
be positive if 


(26, — 1)’p(1 — p) < (26m, — 1)” ap(1 — p,) (22) 


To obtain the condition for @,, satisfying (22) write: 


u=l-—e and mi 
then from (22) 
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developing yu in powers of « 
Ou, > 6, + 4(26, — le + terms in higher powers of «. (23) 


Since ¢ is usually small and as can be seen from (14) the difference 
between @,, and @, increases with ¢ it follows that in most practical 
situations (22) is satisfied when 04, > 6,. Thus when @ = 6) = 04, = 3 
in the majority of cases the matched-pair trial will be shorter than that 
based on unrestricted pairing. Corresponding proofs for Method 2 
are easier to obtain since the common proportionality factor \ simplifies 
the formulae to a considerable extent. 


DISCUSSION 


P. Armitage (1954) in the concluding paragraphs of his paper 
commented on the possibility of using matched pairs in sequential 
trials. Two methods of arranging a sequential plan for matched pairs 
have been discussed here and shown experimentally to give consistent 
results. 

On many occasions matching, even if desirable from the medical 
point of view, may not be practicable. When, however, the flow of 
cases is ample and fairly steady, matching becomes a practical proposi- 
tion. 

In the examples given above, the matched samples were on the 
average about 10% shorter than the series using random pairing from 
the whole population. The differences between strata were, however, 
purposely exaggerated and in practice the difference between total 
sample sizes is likely to be small. In spite of this, matching may be 
advantageous for reasons indicated in the Introduction. 

The formula for @,,, (6) is rather insensitive to moderate errors in 
the estimation of p, and a, so that the method may be used whenever 
reasonably good estimates of these quantities can be obtained. 

If the response in the strata fails to follow the assumed constant 
64 or proportional change pattern the length of the trial will be affected. 
Whether it will be shortened or lengthened will depend on the direction 
of departure from the assumed pattern and on the values of p; and a; 
in the strata affected. In these circumstances, however, the length of 
the trial is of comparatively little importance since we gain valuable 
information on differential response to the experimental treatment. 
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A NOTE ON THE RANK ANALYSIS OF INCOMPLETE BLOCK 
DESIGNS—APPLICATIONS BEYOND THE SCOPE 
OF EXISTING TABLES 


Orro Dykstra, JR.* 
General Foods Corporation, Hoboken, New Jersey, U.S.A. 


1. Introduction 


The use of rank analysis in incomplete block designs using a method 
of paired comparisons is covered in papers by Bradley, either alone 
(1953, 1954a, b, 1955) or in co-authorship with Abelson (1954) and 
Terry (1952). As in these papers, it is assumed that there is an experi- 
ment consisting of ¢ treatments with n repetitions. Furthermore, 
Bradley (1954a) postulates that associated with each of the ¢ treatments, 
denoted by 7, , --- , 7, , there exist parameters, 7; for 7’; , such that 
x; > Oand )-j., x, = 1. The parameters are further defined with the 
probability statement that, if X; generally denotes an observation on 
a sample of 7; , 


(1) P(X; > X;) = + 


in the comparison of 7; with T; , i # j. 
Bradley and Terry give the formulas for the maximum likelihood 
estimates of the =; : 


(2) 
and 

(3) Lp =1, 

and define 

(4) a; = 2n(t — 1) — 


where r;;, = lif X,; > X, andr,;, = 2if X; > X, in the kth repetition 
of the pair (2, 7). 


*I am indebted to Mavis B. Carroll and Clifton C. Sutton for a careful reading of the paper and 
their many constructive comments and criticisms. 
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Tables for the test procedures for small treatment and sample 
sizes are provided in the papers by Bradley and Terry (1952) and Brad- 
ley (1954b). The tables cover ¢ = 3,n = 1 (1) 10; = 4,n = 1 (1) 8; 
and t = 5, n = 1 (1) 5, giving the p; to two decimals. 

We use Equation (2) in the form 


a; / n 
6) Fp) 


Using any method of obtaining initial p; , the resulting first estimates 
are substituted into the right side of (5); the next estimates being 
resubstituted, and so on, until the equalities hold. This is the method 
formerly and still recommended for iterating. The only problem is to 
obtain the initial p; . This iterative procedure must be used even with 
the tables, to obtain more decimals. Some experimentation with the 
procedure will demonstrate the slowness of the convergence to the 
maximum likelihood estimates, the values converging from only one 
side. 

During a seminar at V. P. I. Professor M. G. Kendall suggested 
one method which might give reasonably good first approximations from 
which solutions of the likelihood equation might be found by iteration. 
Letting the number of times 7; > T; in n repetitions of the pair (¢, j) 
be denoted by a,,; , then the first estimates are given by 


(6) = — (t — 2)]. 


_A shortcoming of the formula is that the more the a,; deviate from 
their expected values, the more inaccurate the first estimates given by 
(6) will become. 

The main purpose of this paper is to present a quick and easy 
method of obtaining first estimates of the x, , regardless of the size of 
# and n. In addition, because in certain cases the formula severely 
overestimates some of the first estimates, a means of oquasting for the 
overestimation is also given. 


2. Derivation of the Formula Estimating the x; 


We now develop an alternative formula, which does not have the 
shortcoming of the Kendall formula and which in most experimental 
situations will prove superior to (6) even without the correction to be 
described later. It gives first estimates generally closer to the maxi- 
mum likelihood estimates, thus requiring fewer iterations. 

We assume (and later examine the assumption) that thé p, in (5) 
are not too different from one another. If the p,’s were all equal to 


4 
i 
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p; , say, then using the fact that >>, p, = 1, 


(7) ps = (1 — p)/(t — 1). 
Substituting (7) for the p, in (5) we finally obtain 
a,/n 
which, when solved for p, , yields 


(9) = a,/[n(t — 1)? — aft — 2)}*. 


3. Application of the Formulas for Estimating x, 


The experimental example in Table 1 is given to compare the accu- 
racy of the formulas (6) and (9). In this example ¢ = 8 and n = 40. 


TABLE 1 
Comparison of the Formulas (6) and (9) 


Ps — Se 
t - 
(6)* (9)* (6) (9) 

1 211 208 . 204 — .003 — .007 
2 164 165 .163 001 — .001 
3 147 151 .148 004 001 
4 133 135 .134 002 001 
5 125 137 197 012 002 
6 082 084 .083 002 001 
7 070 056 .072 —.014 002 
8 068 065 .069 — .003 001 


*The values determined by each formula have been adjusted proportionately so that >i ps = 1. 


In this example the a,, did not deviate appreciably from their 
expected values, so that (6) could be applied. The column of differences 
between the two sets of first estimates (p,) and the maximum likelihood 
estimates (#;) indicates the superiority of (9), which would be still 
greater if more decimals were required. 


4. Further Improvement of the Formula p, = a;/{n(t — 1)? — a,(¢ — 2)} 


The advantage in using the equation suggested here for estimating 
the #; is not so great when the approximation based on the assumption 


*These results were also obtained by C. Y. Kramer and M. C. K. Tweedie at the Virginia Poly- 


technic Institute by somewhat different approaches, independent of the present autbor and of each 
other, in unpublished work. 


| ng 
3 
= 
7 
4 
Gk 
ie 
q iy 
re 
iq 
4 
q 
pee 
, 
| 
P 
: 
Agee 
ay 
| 
| 
| 
| 


304 BIOMETRICS, SEPTEMBER 1956 


of equal p,’s is inaccurate (such as occurs when the spread in the various 
x; values is large). ‘This is illustrated by the example in Table 2, 
when ¢ = 4 and n = 144. 


TABLE 2 
Comparison of #; and p; = a;/{n(t — 1)? — ax(t — 2)] 


i a; cp;** pi — — 
1 271 3593 3594 3529 0001 | —.0064 
2 208 2294 2364 2321 0070 0027 
3 193 2064 2121 2083 0057 0019 
4 192 2049 2105 .2067 0056 0018 
Total 1.0000 1.0184 1.0000 .0184 .0000 


*Determined by (9) 
**c adjusts the p: proportionally so that cps = 1. 


It can be seen that each p; in the table overestimates the correspond- 
ing #; , and the amount of overestimation is smallest for the largest p, . 
In general, the spread of the p,’s (j ¥ 7) and thus the amount of over- 
estimation is less for the smallest and largest p,’s than it is for the inter- 
mediate p,’s. Since p; > #; for all 7, >>; p; > 1, so that when the p, 
are reduced proportionately to total one, the largest 2, is severely 
underestimated. 

An amended formula, determined empirically, has been used exten- 
sively and invariably improved the approximation to the x; . In 
most instances indeed the amended formula gave sufficiently accurate 
estimates, with only one iteration necessary to check the convergence 
to the maximum likelihood estimates. We first define | 


(¢ — DR — + aa. 


where 


R=) ai, 


a, =nt— 1) —a,, 
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and 
p; = a,/{n(t — 1)? — a(t — 2)], equation (9). 


Then an improved first estimate (f;) for the maximum likelihood 
estimate of x, is given by: 


(11) B=p—k,. 


5. Application of the Correction for p; = a,/[n(t — 1)? — a,(t — 2)] 


The improvement in the estimates given by (9) when corrected using 
(10) is shown in Table 3. The results of Table 2 are used to illustrate 
the procedure. 


TABLE 3 
Application of (9) and (10); 2 = 4,n = 144 


1 271 161 3594 .0003 3591 3593 |—.0002 
2 208 224 2364 .0068 2296 2294 0002 
3 193 239 2121 .0057 2064 2064 0000 
4 192 240 2105 .0056 2049 2049 0000 

864 864 1.0184 0184 | 1.0000 | 1.0000 0000 


S= = n( 4) = 864 
DR — 


|| 
| 
e & 
or 
& 


It is obvious that the #; are very close to the maximum likelihood 
estimates #;. In this example only one iteration is necessary to obtain 
the #,; using . 

In computing k,; the only term which changes as 7 changes is the 
product a,a,.,; , the other terms being constant for a given experiment. 
When the p, are estimated by (9), the correction need be applied only 
when )., p; is substantially larger than one. _ 
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6. Conclusions 


In summary, we can list the following two advantages for the 
proposed method: 


a. It is quick and easy to apply. Even when the p, are widely 
diverse the two formulas together will give first estimates very 
close to the maximum likelihood estimates. 

b. The first estimates are not sensitive to deviations from the model, 
depending on the sum of the a,, and not the individual a,; , as 
in the Kendall formula. Being insensitive to deviations from the 
model is not a disadvantage, since the maximum likelihood 
estimates of the +, must first be obtained, in order to test the 
appropriateness of the model [Bradley (1954a)]. 


REFERENCES 


Abelson, R. M. and Bradley, R. A. (1954). A2 X 2 factorial with paired comparisons. 
Biometrics, 10, 487. 

Bradley, R. A. and Terry, M. E, (1952‘. Rank analysis of incomplete block designs. 
I. The method of paired comparisons. Biometrika, 39, 324. 

Bradley, R. A. (1953). Some statistical methods in taste testing and quality evalua- 
tion. Biometrics, 9, 22. 

Bradley, R. A. (19542). Incomplete block rank analysis: on the appropriateness of 
the model for a method of paired comparisons. Biometrics, 10, 375. 

Bradley, R. A. (1954b). Rank analysis of incomplete block designs. II. Additional 
tables for the method of paired comparisons. Biometrika, 41, 502. 

Bradley, R. A. (1955). Rank analysis of incomplete block designs. III. Some large- 
sample results on estimation and power for a method of paired comparisons. 
Biometrika, 42, 450. 


+ 
a 


EXTENSION OF MULTIPLE RANGE TESTS TO GROUP 
MEANS WITH UNEQUAL NUMBERS OF REPLICATIONS 


CiypE Youne KRAMER 
Virginia Agricultural Experiment Station of the Virginia Polytechnic Institute 
Blacksburg, Virginia, U.S.A. 

In many fields of research, one is faced with the task of comparing 
the effects of treatments which have been replicated unequally. This 
happens for a number of reasons. In an experiment on animals, some 
may get sick and have to be removed from the experiment. In some 
experiments, the amount of material available for certain treatments 
may not be as much as for other treatments. If the experimenter has 
specified orthogonal contrasts that he is interested in before he runs 
the experiment, one can test the various treatment effects by an F-test 
after the treatment sum of squares has been partitioned into individual 
degrees of freedom for each orthogonal contrast. If the experimenter 
has not specified orthogonal contrasts, one is faced with the problem 
of deciding which treatments are significantly different. 

Several writers, including Duncan, Keuls, Newman, and Tukey, 
have developed multiple range tests to show differences among treat- 
ments that have been ¢eplicated the same number of times and when 
nothing was specified concerning the treatments. Duncan [1] compares 
the above methods and gives citations. This extension to unequal 
numbers of replications will be exemplified with reference to Duncan’s 
“New Multiple Range Test,’”’ but is applicable to any of the above 
writers’ tests; all one has to do is use their tabled ranges. 

In Duncan’s test for an equal number of replications, the differ- 
ence between any two ranked means is significant if the difference 
exceeds a shortest significant range. This shortest significant range is 
designated by R, and is obtained by multiplying the standard error of 
a mean, 8, , by a given value, z,,,, , obtained from a table of significant 
studentized ranges which Duncan has tabled for both the 5% and 1% 
test. In Duncan’s terminology, n, is the degrees of freedom of the error 
mean square and p = 1, 2, --- , ¢is the number of means concerned. 

Consider an experiment with five treatments, A, B, C, D, and E, 
each replicated n times. Suppose on ranking the means from low to 
high one obtains 


Bo £4 Zp Ze Te. 


In order for Z,; — Z to be significant, Z; — Z, must exceed R; = 8,25 .n, ; 
in order for Z; — Z, to be significant, Z; — Z, must exceed R, = 3,%.,,, 
etc. Now if we had only two treatments, B and C, 3 — Z¢ would be 
significant if — exceeded R, = Since = 
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and = s°/n, then 


It is then seen that 


/2 


If 2, , , , are based on , Me , , replications, 
&, = 8'/n,. Now for Zs; — Zc to be significant, it is reasonable that 
— Ze should exceed 


and for Z; — Z, to be significant, Z; — #, should exceed 


+ 


Now e 


1), 
£5 to > 


can be written 

Np + Nc 

which indicates that for grouping means based on unequal numbers, 

a table of factors Rj = sz,,,, , where s is the mean square for error, 

should be set up, in actually making this test in each individual case. 
The above test will be identical to a ¢-test if there are only two means 

since by using Duncan’s test for unequal numbers we’ get 


2, — 20> (2 


(22 — Zc) > 


and by a ¢-test 


—#o>y(L 


if (£g — Zc) is significant. Xx comparing Duncan’s tabled z,,,, with 
we see that = 2. 
As pointed out by a referee, this extension to unequal numbers 
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of replications will be a conservative test. Evaluation of specific 
significance and prediction levels would be extremely difficult and im- 
practicable. If the number of replications differs greatly, there will be 
an increased probability of a significant difference within a subset of 
ranked means classed as homogeneous by this test. 


NUMERICAL EXAMPLE 
(a) Treatment Means Ranked in Order and Number of Replications: 
F D A B Cc E 


z 458 498 6521 528 564 630 

n 3 5 4 3 5 2 
(b) Analysis of Variance: 
Source Di. M.s. F 
Between Treatments 5 9306.17 3.88 (P < .025) 
Error 16 2397 .00 
{c) s = V 2397 = 48.96 
(d) Significant Studentized Ranges for a 5% Duncan Multiple Range 
Test [1]: 

p 2 3 4 5 " 


N = 16 3.00 3.15 3.23 3.30 3.34 
(e) Appropriate Significant Range Factors (Rj = 82,1.) 


p: (2) (3) (4) (5) (6) 
RB: 146.88 154.22 158.14 161.57 163.53 
(f) Test: 
In order for (Zz — Zr) to be significant, 
2)(3 


must exceed Rj = 163.53. 


(Zs — Zr) AOS) = 172 x 1.54 = 264.88 > Ri, 


therefore, — Zp is significant. 
In order for (2g — Zp) to be significant, 


2) 
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must exceed Ri = 161.57. 


(fz — Zp) 245 


5 = 132 X 1.69 = 223.08 > Re, 


therefore, Zz — Zp is significant. 
Completing the test, we have 


— = 109 X 1.68 = 177.67 > Ri, 


therefore, significant, 


(Ep — Zs) see) = 102 X 1.54 = 157.08 > Rf, 


therefore, significant, 


(5) 


(fz — 245 


= 66 X 1.69 = 111.54 


therefore, not significant, 


2(5)(3) 


(te 


= 106 X 1.94 = 205.64 > Rj, 


therefore, significant, 


2018) 


= 66 X 2.24 = 147.84 < Rj 


therefore, not 


2(3)(3) 


= 70 X 1.73 = 121.10 < Ri, 


therefore, not significant. 
Now we may write 


F . D A E 
458 498 521 528 564 630 


and say that any two means not underscored by the same line are 
significantly different, and any two means underscored by the same line 
are not significantly different, in the sense defined. 


REFERENCE 


{1] Duncan, D. B., Multiple Range and Multiple F Tests, Biometrics, 11, 1-42, 1955. 
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SIMPLIFIED (OR ED,;.) CALCULATIONS 


Henry J. Horn . 
Hazleton Laboratories, Falls Church, Virginia, U.S.A. 


In industry a low-cost screening program to evaluate adequately the 
hazard and/or potency of numerous new chemicals, drugs, compounds, 
and combinations thereof is a fundamental requirement. An idea of 
the range of toxicity or potency can be satisfactorily determined from 
a limited number of animals. Four to five animals per group will be 
satisfactory as an accurate evaluation is usually not required in the 
early stages of development. Furthermore, a costly screening program 
is economically unfeasible as from a multitude of new items only one or 
two may reach widespread commercial usefulness. 

In light of the above, this paper presents a practical method for the 
routine screening to estimate the acute LD;) (or ED;.) by whatever 
route of administration is desired. By a unique choice of doses a table 
of LD;. (or ED;o) values with their confidence limits can be prepared; 
hence, the mathematical computation is avoided and all that is required 
is to read the values from the tables. Although only LD,, is spoken of 
throughout the remainder of the paper for convenience of presentation, 
the calculation of the ED;, will be identical. 

The statistical method used is the moving-average interpolation 
method as originally presented by Thompson [1]. In this method, the 
logarithm of the LD;. is determined by simple linear interpolation 
between moving averages of the doses surrounding 50% mortality; 
the antilogarithm will then yield the estimated LD;, . Weil [2] has com- 
puted tables to assist in the computation of the LD;, when four suc- 
cessive doses (D, , D, , D; , Ds) in geometric progression are used and 
the same number of animals (n) is used at each dosage level: 
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(1) Log LDs. = log D, + d(f + 1) 
where: D, = lowest dose of the 4 doses used in the calculation 

d = logarithmic increment between logarithms of the 
doses or logarithm of the multiplier relating each 
successive dose to its immediately smaller dose. 

f = tabled value given by Weil [2] for the combination 
of r, , 2 , 73 , and r, (actual number of animals 
responding, for example, dead, at doses D, , D; , 
D, and D, , respectively) determined experiment- 
ally. 


Further simplification of Equation (1) can be accomplished as the 
logarithmic increment is constant between doses; therefore, the one d 
may be dropped when D, is substituted for D, . Furthermore, the 
equation may be written in terms of the LD, rather than the logarithm 
of the LD,.. This can be accomplished by rewriting, as follows: 


(2) LD = antilog [log D, + d X f] 
The value of f will vary from 0 to 1; 
when f = 0 then the LD,. = D, and 
when f = 1 then LD, = D;. 


The 95 per cent confidence limits may be calculated by the following 
surmula: 


(3) Coafidence limit = antilog (log LD,. + 2 X d Xa) 


‘fables of values of o, are given by Weil [2]. If the value of a, is zero, 
whieh occurs only when r, , r2 , 7; , and 7, are zero mortality, zero 
mortality, total mortality, and total mortality, respectively, the con- 
fidence limits cannot be calculated mathematically but will lie within 
the boundaries specified by those successive doses the lower of which 
produced no mortality while the higher produced total mortality. 
The main purpose of this paper is to present a further simplification 
of the moving average interpolation method. This is accomplished by 
using a dosage increment factor which is an integer root of n; thus the 
logarithmic increment between the logarithm of doses will be a simple 
fraction. Examples are 0/10, ~/10, etc., and the logarithms are 
1/2, 1/3, etc., respectively. By this choice, it will be found that the 
series of doses will repeat themselves at successively higher multiples 
of 10; for example, the dosage series for a multiplier of 1/10 may be 
expressed as ee X 10‘ where ¢ = 0, +1, +2, +3, etc. This means 


that the mantissae of the logarithms will be found to repeat themselves 
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at intervals of every other one, two, etc., depending on whether V 10, 


10, etc., is used, respectively. 


The following series of practical doses are proposed for the estima- 


tion of the LD, : 


Multiplier = 710 
Log increment = 1/2 


. $1.00 
Dosage series xX 10 


Multiplier = ~/10 
Log increment = 1/3 


1.00 
Dosage series {2.157 X 10° 


4.64 
Dosage level, Dosage level, 
mg/kg mg/kg 
10.0 10.0 
31.6 21.5 
100 46.4 
316 100 
1000 215 
3160 464 
1000 
2150 
4640 


It will be noted that doses range from a low of 10 mg/kg to a high of 3 
to 5 g/kg; this is adequate for screening purposes. If a compound has 
an LD, below 10 mg/kg it will be known to be extremely toxic and, 
conversely, if the LD, is greater than 3 to 5 g/kg, it can be assumed 
to be relatively nontoxic. 

As there is a cyclic repetition of dosage increment levels, there will 
also be a cyclic progression of corresponding LD,. values. Hence, if a 
table of the LD,, is prepared for one cycle, say between 1 and 10, then 
the values for any other cycle will be 10‘ times the tabulated values 
(where ¢ = 0, +1, +2, +3, etc.). The LD, values must be tabulated 
according to the doses (D, and D,) that straddle the LD, ; that is, 
with a dosage series a < 10‘, two columns of LD,o values are 
‘presented, the first column covering those which lie between 1.00 X 10° 
and 3.16 X 10‘, while the second column covers those which lie between 
3.16 X 10‘ and 10.0 X 10‘. As the multiplier becomes a higher root 
of 10, the number of doses per cycle increases; hence, an increasing 
number of columns of LD,. values must be listed (equal in number to 
the power of the root). 
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Four tables of LD,;. values and their 95 per cent confidence limits 
have been prepared: Table 1 for n = 4 (4 animals at each dosage level) 
and a multiplier of V/10 between successive doses; Table 2 for n = 4 
and a multiplier of ~/10; Table 3 for n = 5 and a multiplier of V10; 
and Table 4 for n = 5 and a multiplier of ~/10. These tables have 
been computed using the values of f and a, as given by Weil [2]. 

In practical application, the series of dosages specified in the pre- 
vious lists will be administered. Selection of the number of animals 
per dose (i.e., 4 or 5) and the particular series of doses (i.e., using a 
multiplier of ‘V10 or V/10) will be determined usually by economics. 
If the range of toxicity is to be determined in larger animals such as 
rabbits, cats, etc., four animals per dosage level and the series correspond- 
ing to the 10 will be more economical. However, if smaller animals 
such as mice or rats are used, the use of five animals per dosage level 
and a series of dosages corresponding to W10 will usually be found 
feasible. 

An experienced technician can also effect further improvement in 
economy by selectively administering only those dosage levels which 
are necessary for the determination of the LD, . For example, if the 
LD,. is below 10 mg/kg the administration of only the low dosage 
levels will be adequate; conversely, if the LDzg. is above 5 g/kg, the 
administration of the higher dosage levels only will be adequate. Fur- 
thermore, it is often possible to make a rough guess as to the LD,. and 
administer only those dosage levels in that vicinity. If there is no 
knowledge of the LD;, , then experience has shown that administration 
of the middle doses will often ciarify the situation. -For example, if all 
the animals at say 215 mg/kg die within 5 to 10 minutes, then ony 
the lower dosage levels need be given. 

Occasionally the mortality results of the-series of dosages will yield 
no mortality for the lower doses followed by an immediate transition to 
total mortality for the successive higher doses. This means that for 
the particular compound under study, the regression line is relatively 
steep. Under these circumstances, it will be impossible to determine 
confidence limits although it will be quite obvious that these values 
lie between those doses at the transition point. If, or when, a more 
accurate LD,» need be determined, the range of doses for the new study. 
will be apparent. 

In addition to the convenience of determining the LD; , other 
advantages arise from the use of this scheme of doses. The preparation 
of solutions containing the test material will be simplified since 1 to 10 
and 1 to 100 dilutions will yield the correct dosages at the lower levels 
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when the same quantity per kg. is administered; that is, the factor for 
the quantity of material to be administered per unit of body weight 
will remain the same in each dosage cycle. It will also tend to maintain 
the quantity administered to all the animals within closer limits. 

As only four or five animals are used per dosage level, the selection 
of animals for each group in the study is most important. Because 
individual variations may be marked, poor selection will produce an 
irregular progression in the mortality results. In order to keep this 
irregular progression to a minimum, a stratified random sampling of 
animals is recommended. This is obtained, for example, by placing in 
each dosage group equal numbers of individuals having each of the 
following characteristics: heaviest weight, lightest weight, average 
weight, most agile, most lazy, etc. Randomization is used to allocate 
the animals from each category to all the dosage groups. Animals of 
questionable health must be excluded. 

The following example will serve to illustrate the use of the simplified 
technique presented here: It was desired to determine the rat oral LD;, 
of an unknown compound on which nothing was known relative to 
toxicity. Five rats per group were used with the dosage progression 


4.64 
was administered first and after one-half hour, signs of ataxia were 
observed. It was therefore felt that administration of the next higher 
dose and all the lower doses would result in bracketing the LD5 . 
Results at the end of seven days were: 


1.00 
past X 10° as described above. The middle dosage of 215 mg/kg 


Dose Mortality 
mg/kg per group 
10.0 0/5 
21.5 0/5 
46.4 1/5 
100 2/5 
215 4/5 

464 5/5 


Inspection indicates that either the series of doses with corresponding 
mortality of 0, 1, 2, 4, or that with 1, 2, 4, 5 should be analyzed. A 
search through Table 4 shows that the series 1, 2, 4, 5 is listed and that 
the estimated LD,, is 110 mg/kg with confidence limits of 55.0 to 220 
mg/kg. 
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SUMMARY 


The moving average interpolation method has been summarized. 
A far more convenient and rapid use of this method may be achieved 
by choosing a dosage increment factor which is an integer root of 10; 
e.£g., V10, W10. This produces a repetition of the dosage levels at 
successive multiples of 10 and simultaneously allows for the prepara- 
tion of tables giving estimated LD,. or EDs. and their confidence 
limits directly from the observations. Four tables have been prepared 
using parameters which experience has shown to be most useful. It 
must be emphasized that this method is not applicable for an accurate 
determination of the LD; or ED,. , but, rather, for estimating the 
range of toxicity and/or potency as in screening programs. For this, 
the economic advantage far overshadows any disadvantage due to 
inaccuracy. 


REFERENCES 
{1] Thompson, W. R. Use of moving averages and interpolation to estimate median 
effective dose. Bact. Rev., 11, 115-145, 1947. 


[2] Weil, C. S. Tables for convenient calculation of median-effective dose (LDw or 
ED, ) and instructions in their use. Biometrics, 8, 249-262, 1952. 
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LDw (OR EDs) CALCULATIONS 


TABLE 1 


= x 
3 
~ 
oD 
H 
= w 
- 
= 
AEs 
gx 
= 
que. € & 
= 
Q 


Confidence 


LDso 
10.0 


limits 


LDso 
3.16 
2.37 
1.78 


— 19.4 


5.14 


6.15 
4.22 


4.22 — 13.3 


7.50 


1 
4.01 


0.750 — 2.37 


84 
0.788 — 


noe 
SAM 


an 
se 
wo 
N 


OTN 


oon N 


1.40 — 15.4 


4.64 
3.16 


2.96 


0.412 — 2.43 


0.338 — 


0.442 — 4.87 


1.47 
1.00 
1.00 


ooo 


0.621 — 16.1 
0.788 — 


3.16 
1.78 


2.49 — 12.7 
1.22 — 81.8 


0.911 — 34.7 


10.0 


5.62 
3.16 
5.62 
3.16 


0.288 — 11.0 


1.78 


0.621 — 16.1 


0.725 — 43.6 


0.229 — 13.8 


1.78 


0.386 — 25.9 


0.221 — 45.2 


3.16 


0.0700— 14.3 


1.00 


7 
| 
he 
| 
0 0 2 4 1.63 — = 
0 0 3 4 1.33 — 
1 1 4 3.16 1.40 — 10.0 4.43 — 22.6 
1 2 4 2.37 | 7.50 3.11 — 18.1 pike 
1 3 4 1.78 5.62 2.49 — 12.7 ieee 
1 4 4 1.33 4.22 2.37 — 7.50 ene 
2 4 3 1.30 — 7.67 
3 3 3 a 3.16 1.07 — 9.36 : es 
1 0 3 3 10.0 1.96 — 50.9 ; 
1 1 2 3 3.16 0.386 — 25.9 —_ ae 
1 1 3 3 
1 1 4 3 1.00 0.196 — 5.09 tag 
1 2 2 3 
1 2 3 3 1.00 0.122 — 8.18 
2 0 3 3 3.16 0.122 — 82.1 10.0 0.385 —260. 
2 0 4 3 
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TABLE 1 (continued) 


7 Di = 0.316 Di= 1.00 
| Dz: = 1.00 Da = 3.16 
or *™ pone *™ 
T1,73,72,76 Ds = 10.0 Ds = 31.6 
Confidence Confidence 
LDso limits LDso limits 
2 1 2 3 3.16 0.0472—212. 10.0 0.149 —669. 
2 1 3 3 1.00 0.0149— 66.9 3.16 0.0472—212. 
2 2 2 3 1.00 0.0100—100. 3.16 0.0316—316. 
0 0 4 2 3.16 0.837 — 11.9 10.0 2.65 — 37.8 
0 1 3 2 3.16 0.386 — 25.9 10.0 1.22 — 81.8 
0 1 4 2 1.78 0.471 — 6.72 5.62 1.49 — 21.2 
0 2 2 2 3.16 0.316 — 31.6 10.0 1.00 —100. 
: 0 2 3 2 1.78 0.2971 — 11.7 5.62 0.858 — 36.9 
oa 0 2 4 2 1.00 0.265 — 3.78 3.16 0.837 — 11.9 
“a 0 3 3 2 1.00 0.196 — 5.09 3.16 0.621 — 16.1 
7 1 0 4 2 3.16 0.221 — 45.2 10.0 0.700 —143. 
1 1 3 2 8.16 0.0472—212. 10.0 0.149 —669. 
1 1 4 2 1.00  0.0385— 26.0 3.16 0.122 — 82.1 
1 2 2 2 3.16 0.0316—316. 10.0 0.100—1000. 
1 2 3 2 1.00 0.0149— 66.9 3.16 0.0472—212. 
: 0 2 3 1 3.16 0.0472—212. 10.0 0.149 —669. 
bi 0 2 4 1 1.00 0.0700— 14.3 3.16 0.221 — 45.2 
% 0 3 3 1 1.00 0.0385— 26.0 3.16 0.122 — 82.1 
~ 0 1 4 1 3.16 0.122 — 82.1 10.0 0.385 —260. 
ace TABLE 2 
5 Values of estimated LD,, and confidence limits for geometric dosage 
4 1.00 


series (2.15) X 10' where = 0, +1, +2, +3, ete. (multiplier = ¥/10). 
4.64 
Four animals per dosage level. 


Di = 0.464 Di = 1.00 Di = 2.15 
r1,72,73,74 Dz = 1.00 Dz = 2.15 Di = 4.64 
or preset pesal*™ nena *™ 
Ds = 4.64 De = 10.0 Ds = 21.5 
Confidence Confidence Confidence 
LDso limits LDso _—limits LDso limits 
: 0 0 2 4 2.15 1.388 —3.36 4.64 2.98— 7.23 10.0 6.42 — 15.6 
0 0 3 4 1.78 1.21 — 2.61 3.83 2.61 — 5.62 8.25 5.62 — 12.1 
0 0 4 4 1.47 _ 3.16 _ 6.81 — 
0 1 1 4 2.15 1.25 — 3.71 4.64 2.70 — 7.99 10.0 5.81 — 17.2 
0 1 2 4 1.78 0.989 — 3.20 3.83 2.13 — 6.89 8.25 4.59 — 14.8 
0 1 3 4 1.47 0.853 — 2.53 3.16 1.84 — 5.44 6.81 3.96 — 11.7 
0 1 4 4 1.21 0.825 — 1.78 2.61 1.78 — 3.83 5.62 3.83 — 8.25 
: 0 2 2 4 1.47 0.784 — 2.75 3.16 1.69 — 5.92 6.81 3.64 — 12.7 
. 0 2 3 4 1.21 0.674 — 2.18 2.61 1.45— 4.69 5.62 3.13 — 10.1 
| 0 2 4 4 1.00 0.642 — 1.56 2.15 1.38 — 3.96 4.644 2.98 — 7.23 
| 0 3 3 4 1.00 0.581 — 1.72 2.15 1.25— 3.71 4.64 2.70 — 7.99 
{ 1 0 2 4 2.15 1.19 — 3.89 4.64 2.57 — 8.38 10.0 5.54 — 18.1 
1 0 3 4 1.67 0.973 — 2.86 3.59 2.10— 6.16 7.74 4.51 — 13.3 
1 0 4 4 1.29 0.918 — 1.82 2.78 1.98 — 3.91 5.99 4.26— 8.43 
1 1 1 4 2.15 1.04 — 4.44 4.64 2.25— 9.57 10.0 4.85 — 20.6 
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Confidence 


limits 


5 
1.70 — 58.9 


4.12 — 24.3 
2.81 — 16.5 
4.64 0.529— 40.7 
4.85 — 20.6 
3.93 — 25.5 
3.38 — 29.6 
2.46 — 40.6 
1.74 — 26.7 


Di = 2.1 
Di = 4.64 
Ds = 10.0 
De = 21.5 
LDso 
7.74 3.48 — 17.2 
5.99 2.69 — 13.3 
4.64 2.25— 9.57 
5.99 2.43 — 14.8 
4.64 1.82 — 11.8 
10.0 
10.0 
4.64 1.14 — 18.8 
4.64 1.00 — 21.5 
10.0 
4.64 0.282— 76.5 
10.0 
10.0 
7.74 3.48 — 17.2 
7.74 3.14 — 19.1 
5.99 2.69 — 13.3 
10.0 
6.81 3.96 — 11.7 
10.0 
6.81 2.02 — 22.9 
4.64 1.57 — 13.7 
6.81 
4.64 1.14 — 18.8 
4.64 0.789— 27.3 
10.0 0.607—165. 


8.00 
6.19 
5.44 
8.75 


4.64 0.529— 40.7 
2.15 0.366— 12.7 
4.64 0.282— 76.5 
2.15 0.131— 35.5 
2.15 0.100— 46.4 


4.64 1.91 — 11.3 


8.00 
6.19 
4.44 


6.86 


2.15 0.846— 5.48 


4.64 1.91 — 11.3 
3.16 1.30 — 


9.57 
8.86 


limits 


Confidence 
3.59 1.61 — 


1.57 — 13.7 


3.16 0.903— 11.1 


2.25 — 


2.15 0.888— 5.23 
2.15 0.728— 6.38 


2.15 0.531— 8.75 
2.15 0.464— 10.0 
3.16 0.807— 12.4 


4.64 0.789— 27.3 


2.15 0.246— 18.9 
4.64 0.529— 40.7 
2.15 0.131— 35.5 
4.64 1.57 — 13.7 
4.64 1.14 — 18.8 
3.16 0.940— 10.6 


4.64 1.82 — 11.8 
3.16 1.84 — 


2.15 0.531— 


3.59 1.46 — 
2.78 1.25 — 


2.78 1.25 — 
3.59 1.61 — 


2.15 1.04 — 


2.78 1.13 — 
4.64 
4.64 


LDso 


TABLE 2 (continued) 


Confidence 


limits 


1.67 0.750 — 3.71 
0.524 — 3.18 
0.246 — 4.06 
0.215 — 4.64 
0.728 — 6.38 


0.485 — 2.06 


2.15 0.888 — 5.23 
1.47 0.605 — 3.56 
1.00 0.412 — 2.43 
2.15 0.728 — 6.38 
1.47 0.419 — 5.14 
2.15 0.366 —12.7 

1.00 0.114 — 8.77 
2.15 0.246 —18.9 

1.00 0.0607—16.5 

2.15 1.04 — 4.44 
2.15 0.846 — 5.48 
1.67 0.676 — 4.11 
2.15 0.031 — 8.75 
1.47 0.436 — 4.94 
1.00 0.338 — 2.96 
1.47 0.375 — 5.75 
1.00 0.246 — 4.06 
2.15 0.246 —18.9 

1.00 0.170 — 5.89 
2.15 


Di = 0.464 
Dz = 1.00 
Ds = 2.15 
Ds = 4.64 
LDso 
1.00 
1.00 
2.15 


4 
4 
4 
3 
3 


or 


T1,72,73,74 


(OR EDs) CALCULATIONS 


0.131 —35.5 


N 


0.215—100. 
4.12 — 24.3 
2.46 — 40.6 
2.81 — 16.5 
2.15 — 46.4 


0.282— 76.5 


4.64 
4.64 


6.81 


10.0 
10.0 
10.0 


1.14 — 18.8 
3.16 1.30 — 7.67 
1.00 — 21. 


4.64 


4.64 


2.15 0.531 — 8.75 
1.47 0.605 — 3.56 


1.95 — 23.9 
1.91 — 11.3 


| 
= 


1.57 — 13.7 
0 .607—165. 
4.64 0.529— 40.7 
0 .464—215. 
0.607—165. 


10.0 
4.64 0.282— 76.5 


10.0 
4.64 0.789— 27.3 


4.64 0.529— 40.7 


6.81 
10.0 


4.64 
4.64 


10.0 


5.23 


3.16 0.903— 11.1 
2.15 0.728— 6.38 
4.64 0.789— 27.3 
4.64 0.282— 76.5 
0.246— 18.9 
4.64 0.215—100. 
2.15 0.131— 35. 
4.64 0.282— 76.5 
0.366— 12.7 
2.15 0.246— 18.9 
4.64 0.529— 40.7 


2.15 0.888— 


2.15 
2.15 


1.47 0.419 — 5.14 
1.00 0.412 — 2.43 
1.00 0.338 — 2.96 
2.15 0.366 —12.7 
2.15 0.131 —35.5 
1.00 0.114 — 8.77 
2.15 0.100 —46.4 
1.00 0.0607—16.5 
2.15 0.131 —35.5 
1.00 0.170 — 5.89 
1.00 0.114 — 8.77 
2.15 0.246 —18.9 


1 


3 


|| 
Di = 1.00 
D2 = 2.154 x 108 
Di = 4.64 
Di = 10.0 
2 0 2 gf | 
2 0 3 7 
2 0 4 
21 2 
2 1 3 
2 2 2 
30 2 
3 0 3 
3 1 1 
0 0 38 
0 0 4 
1 2 
o 1 4 
& 
1 0 3 
1 0 4 a 
1 1 4 
1 2 2 
#1 2 8 
1.00 0.0607—16.5 ah 
22 2 8 1.00 0.0464—21.5 
0 0 4 2 2.15 0.88 —523 
0 2 2 2 2.15 0.464 10.0 . carat 
2 3 2 
0 3 3 2 
1 0 4 2 
12 2 2 
.2 «4 
2 4 4 
} 
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TABLE 3 
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0, +1, +2, +3, etc. (multiplier 


\ 10° where ¢ 


Five animals per dosage level. 


1.00 
3.16 


xf of estimated LD,, and confidence limits for geometric dosage 


limits 


5.07 —15.7 


Confidence 


Di = 1.00 
D: = 3.16 
D: = 10.0 
Ds = 31.6 


LDso 
8.91 


X 10 
Confidence 
limits 
1.60 — 4.95 


Di = 0.316 
D: = 1.00 
Ds = 3.16 
Ds = 10.0 


LDso 


or 
r1,73,72,T4 


ri ,72,73,74 


2.82 
2.24 


3.42 —14.7 
2.93 —10.8 


4.47 —11.2 
4.30 —18.5 


7.08 
5.62 
8.91 
7.08 
5.62 
4.47 
7.08 
5.62 
4.47 
3.55 
4.47 
3.55 
8.66 
6.49 
4.87 


8.66 


1.41 — 3.55 


1.36 — 5.84 
1.08 — 4.64 


0.927— 3.41 


2.82 
2.24 
1.78 
1.41 
2.24 
1.78 
1.41 
1.12 
1.41 
1.12 
2.74 


2.82 — 7.08 
3.19 —15.7 


2.53 —12.5 


0.891— 2.24 


oo 


1.71 — 7.35 


2.16 — 9.25 
2.02 — 6.24 
2.01 — 9.92 
4.26 —17.6 
3.51 —12.0 
3.40 — 6.98 
3.48 —21.6 
2.55 —16.5 


2.00 —11.9 


6.49 
4.87 
3.65 


6.49 


1.01 — 4.97 


0.801— 3.95 
0.638— 1.97 
0.636— 3.14 
0.542— 2.32 
1.35 — 5.56 
1.07 — 2.21 
1.10 — 6.82 
0 .632— 3.75 


0.682— 2.93 


1.54 
2.74 


1.54 


1.70 — 7.85 
2.34 —18.0 
1.60 —14.1 


1.29 —10.3 
1.20 —11.2 
3.18 —21.4 
+2.29 —13.8 
1.75 — 8.39 
2.43 —28.1 


0.740— 5.70 


2.05 


1.53 —20.7 
1.00 —14.6 
1.37 —23.0 


0.819—17.9 


Nn 


4.03 —18.6 
4.53 — 9.31 


0.259— 5.67 


1.21 


2.67 —15.8 
2.63 — 9.01 
2.83 —26.5 


2.25 —18.7 


3.06 —24 


£8853 


|| 
0). 
: 
0 0 4 
0 0 5 5 1.78 
| 
1 4 5 
1 5 5 
0 2 2 5 Pe 
0 2 3 5 Pe 
0 2 4 5 
0 2 5 5 
4 0 3 3 5 
0 3 4 5 
“4 1 0 3 5 
= 1 0 4 5 2.05 
1 0 5 5 
1 1 5 | 
1 1 3 5 2.05 
1 1 4 5 a 
1 1 5 1.15 0.537— 2.48 
1 2 2 5 
; 1 2 3 5 1.54 0.534— 4.44 4.87 
‘ 1 2 4 5 1.15  0.408— 3.27 3.65 
7 1 3 3 5 1.15 0.378— 3.53 3.65 
2 0 3 5 2.61 8.25 
tT; 2 0 4 5 1.78 0 5.62 
a 2 0 5 5 1.21 0 3.83 
2 1 2 5 2.61 8.25 
: 1 3 5 1.78 5.62 
2 1 4 5 1.21 TT | 3.83 P| 
2 5 1.78  0.434— 7.28 5.62 
2 2 3 5 3330 
0 0 4 4 2.74 1.27 — 5.88 8.66 P| 
; 0 0 5 4 2.05 1.43 — 2.94 6.49 
0 1 3 4 2.74 8.66 
0 1 4 4 2.05 6.49 
0 2 2 4 2.74 8.66 
* 0 2 3 4 2 6.49 
4 0 2 4 4 0 3 4.87 1.91 —12.4 
0 2 5 4 0 2 3.65 1.80 — 7.42 
0 3 3 4 0 4 4.87 1.76 —13.5 
3 4 4 0 2 3.65 1.47 9.10 
1 0 4 4 0 7 8.25 3.01 —22.6 
1 0 5 4 1 3 5.62 3.27 — 9.68 
1 1 3 4 2.61 0.658—10.4 8.25 2.08 —32.7 
2 1 1 4 4 1.78 0.528— 5.98 5.62 1.67 —18.9 
4 1 1 5 4 1.21 0.442— 3.32 3.83 1.40 —10.5 
a 1 2 2 4 2.61 0.594—11.5 8.25 1.38 —36.3 


LDs (OR EDs) CALCULATIONS 


TL, 72,73,76 
or 
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TABLE 3 (continued) 


Di = 0.316 

Dz = 1.00 

| *™ 

Ds = 10.0 

Confidence 

LDso limits 
4 1.78  0.423— 7.48 
4 1.21 0.305— 4.80 
4 1.21 0.276— 5.33 
4 2.37 0.539—10.4 
4 1.33 0.446— 3.99 
4 2.37 0.307—18.3 
4 1.33 0.187— 9.49 
4 2.37 0.262—21.4 
4 1.38 0.137—13.0 
3 2.61 1.19 — 5.71 
3 2.61 0.684— 9.95 
3 1.78 4.37 
3 2.61 0.558—12.2 
3 1.78  0.484— 6.53 
3 1.21 0.467— 3.14 
3 1.78 0.434— 7.28 
3 1.21 0.356— 4.12 
3 2.37 0.793— 7.10 
3 2.37 0.333—16.9 
3 1.33 0.303— 5.87 
3 2.37 0.244—23.1 
3 1.33 0.172—10.3 
3 1.33 0.148—12.1 

TABLE 4 


Di = 1.00 
D2: = 3.16 
Ds = 10.0 
Ds = 31.6 
Confidence 
limits 
1.34 —23.6 
0.966—15 .2 
83 0.871—16.8 


X 10¢ 


50 0 .830—67 .8 
22 0.4 1.0 
25 3.77 —18.1 
25 2.16 —31.5 
62 2.29 —13.8 
1.76 —38.6 
62 1.53 —20.7 
83 1.48 — 9.94 
62 1.37 —23.0 
83 1.13 —13.0 
50 2.51 —22.4 


22 0.958—18 .6 
-50 0.771—73.0 
-22 0 .545—32.6 
-22 0.467—38.1 


uo 


Values of estimated LD,, and confidence limits for geometric dosage 


1.00 
2.15 
4.64 


series 


X 10° where ¢ = 0, +1, 


£2, 


10). Five animals per dosage level. 


Tl, 72,73,76 
or 
T1,73,7T2,74 


KKK OOS 
Dar wrwarwndar wd & 


Di = 0.464 
D: = 1.00 
D, = 2.15 * 
Ds = 4.64 
Confidence 
LDso limits 
2.00 1.37 — 2.91 
1.71 1.26 — 2.33 
1.47 
2.00 1.23 — 3.24 
1.71 1.05 — 2.78 
1.47 0.951 — 2.27 
1.26 0.926 — 1.71 
1.71 1.01 — 2.91 
1.47 0.862 — 2.50 
1.26 0.775 — 2.05 
1.08 0.741 — 1.57 
1.26 0.740 — 2.14 
1.08 0.665 — 1.75 
1.96 1.22 — 3.14 
1.62 1.07 — 2.43 
1.33 1.05 — 1.70 
1.96 1.06 — 3.60 


Di = 1.00 
D: = 2.15 
Ds = 10.0 
Confidence 
LDso limits 
4.30 2.95 — .26 
3.69 2.71 — 5.01 
3.16 
4.30 2.65— 6.98 
3.69 2.27— 5.99 
3.16 2.05 — 4.88 
2.71 2.00 — 3.69 
3.69 2.17 — 6.28 
3.16 1.86 — 5.38 
2.71 1.67 — 4.41 
2.33 1.60 3.39 
2.71 1.59 — 4.62 
2.33 1.43 — 3.78 
4.22 2.63 — 6.76 
3.48 2.31 — 5.24 
2.87 2.26 3.65 
4.22 2.29 — 7.75 


+3, etc. (multiplier = 


Di = 2.15 
Di = 4.64 
| 
De = 21.5 
Confidence 
LDso limits 
9.26 6.36 — 13.5 
7.94 5.84 — 10.8 
6.81 
9.26 5.70 — 15.0 
7.94 4.89 — 12.9 
6.81 4.41 — 10.5 
5.84 4.30 — 7.94 
7.94 4.67 — 13.5 
6.81 4.00 — 11.6 
5.84 3.60 — 9.50 
5.01 3.44 — 7.30 
5.84 3.43 — 9.95 
5.01 3.08 — 8.14 
9.09 5.66 — 14.6 
7.50 4.98 — 11.3 
6.19 4.87 — 7.87 
9.09 4.94 — 16.7 


por 
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A SIMPLE METHOD FOR FITTING AN ASYMPTOTIC 
REGRESSION CURVE 
H. D. Patrerson 
Rothamsted Experimental Station, Herpenden, Herts., England 
Introduction 


The asymptotic regression equation 
y=a-— Bp where 0<p<1l 


is frequently used in agricultural research to express the relationship 
between the yields of a crop, y, and the amounts of a given fertilizer, 
x, applied to that crop. This equation is slightly different to that used 
by Stevens (1951) in that a negative sign is given to the term Bp’. 
This change allows y to increase with positive values of 8. 

Stevens (1951) has described an efficient method for determining 
the parameters a, 8 and p. He provided tables for 5, 6 and 7 equally 
spaced ordinates which considerably reduce the arithmetic labour of 
his method and which have the valuable feature that they enable the 
variances and covariances of the estimates of a, 8 and p to be directly 
determined. These tables have been extended by 8. Lipton on the 
electronic computer at Rothamsted to cover the range of n = 3 to 
nm = 12. Since Stevens’ method is one of successive approximation, 
it is admirably suited for use on the electronic computer. 

Pimentel Gomes (1953) has shown that, with equally spaced ordi- 
nates, efficient estimates of p can be obtained by solving equations of 
the type 


Joyo + + = 0 (1) 


where the J are complicated polynomials in r, the estimates of p. He 
also provided tables of the numerical values of the polynomials for 
the cases n = 4 and n = 5. In practice r can be obtained quickly 
using these tables by a process of trial and error. 

The present paper deals with simple methods of estimation, for the 
cases n = 4, 5, 6 and 7, which require the use neither of complicated 
machinery nor extensive tables. The methods suggested are not fully 
efficient (except in special circumstances) but the loss of precision is, 
in most cases arising in practice, small. 
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Only the estimation of p will be considered in detail as the estimates 
of a and @ are given by the linear regression of y on r*, where r is the 
estimate of p. 

For convenience the ordinates will be taken to be 0, 1, 2 --- n — 1. 
The estimates r will, therefore, normally require to be converted to r'”* 
where s is a single increment of fertilizer in standard units. 

_The method of deriving the estimates considered in the present paper 
depends on the result that the efficient estimate of p can always be 
expressed as the ratio of two contrasts between the y with coefficients 
given by functions of p. It can be expected that the ratio obtained by 
using coefficients appropriate to some value of p, equal to po say, will 
be most efficient for a range of values of p around p, . In fact pp can 
be selected so that the whole useful range of p is covered with reasonable 


efficiency by a single set of coefficients for each of the four cases con- 
sidered. 


Estimation of p with four equally spaced levels of fertilizer 


With three levels of fertilizer full information on p is provided by 
the estimate 


Yi — Yo 
This equation can conveniently be rewritten 
Yo (L+ny+ry = 0. (2) 


When a fourth level of fertilizer is included r can still be determined 
from Equation (2), or from ; 


ys (li = 0, (3) 


or from a combination of Equations (2) and (3). Thus the sum of 
Equations (2) and (3) leads to the estimate 


(4) 
Y2 — Yo 


which has frequently been recommended. 

It is clear, however, that the relative amounts of information on p 
provided by Equations (2) and (3) depend on the true value of p. Thus, 
if p is small the range of values from y, to y2 will be greater than that 
from y, to y; and Equation (2) will provide more information than 
Equation (3). If, on the other hand, p is nearly 1 so that the responses 
to each increment of fertilizer are almost equal, Equations (2) and (3) 
will provide approximately equal amounts of information. 
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_ The equations can, therefore, be combined with the relative weights 
u: 1 where yp is to be determined: 


Ys + (u— 7 — — (ut ur + ury = 0. (5) 
The estimate of r given by this equation is 


+ — — BY. (6) 
Yo + — — 
In order to utilise all the information provided by Equations (2) 
and (3) u should be chosen so that the variance of r is minimised. The 
required expression for yu is 


(7) 
+ 3p + 4p +1 

It is interesting to note that with these values of u the estimates 
(6) are fully efficient, i.e. the same as those obtained by the methods of 
Stevens (1951) and Pimentel Gomes (1953). A proof of this statement 
will be given when higher values of n are considered. 

The values of » range from » = 2 when p = 0 tow = 1 when p = 1. 
When » = 1 equation (6) reverts to equation (4). The estimate 
(Ys — ¥:)/(Y2 — Yo) is therefore likely to be suitable only for high values 
of p, with efficiency decreasing with p. Similarly (y; + y. — 2y:)/(y2 + 
Y, — 2yo) will be most efficient when p is small. 

It seems reasonable to consider the estimate (6) when an intermediate 
value of pis used. A value of u = 1.25, corresponding to p approximately 
equal to 0.26, has been chosen. The estimate 


4ys + — Sy: 

= 8 

has an efficiency of over 95% over the whole practical range of values 
of p, as shown in Table 1. 


Estimation of p with five, six, or seven equally spaced levels 
For n = 5, 6 or 7 p can be estimated in a similar manner from the 
ratio of two simple contrasts between the values of y. Additional . 
equations such as 
¥— (1 + ry. = 0 (9) 
are combined with Equations (2) and (3) using weights 1, y, , ue ete. 
The general form of the estimates obtained in this way is 
Yn-2 + — 1)yn-3 + (us = M1) Yn-4 — Mn—4)Y1 — Mn-3Yo 
(10) 
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TABLE 1 
Percentage Efficiencies of Proposed Estimates of p 


Equa- 
~~ (8) (17) (18) (19) 
p n=4 n=5 n=6 n=7 

0.0 89.3 77.4 65.2 60.0 
0.1 95.4 88.5 79.2 74.5 
0.2 98.6 95.7 90.2 86.8 
0.3 99.8 99.2 97.0 95.0 
0.4 99.9 99.9 99.6 98.1 
0.5 99.4 98.9 98.8 96.8 
0.6 98.5 97.0 96.1 92.4 
0.7 97.4 94.7 92.5 86.9 
0.8 96.3 92.3 88.8 81.3 
0.9 95.3 90.0 85.2 76.3 
1.0 94.2 87.8 81.9 72.0 


Values of the » can be chosen so that these estimates are the same as 
those given by Stevens. For the purpose of proving this statement it is 
convenient to rewrite Equation (10) as 


ALY: 0, (11) 
where the \, are subject to the restrictions ‘ 
>> a. = 0, = 0, and = 1. 


The first two restrictions arise because there are n values of \, but 
only n — 3 values of u. 
If the \, are chosen so that 


= Fi, + Fur? + (12) 

Equation (11) and the three restrictions become 
Fe Dy Dy +F, =0, (13) 
nF., +F,, Da" =0, (14) 
j+% +F.> =0, (19 
and F., + + F,, = 1. (16) 


Stevens has shown that his estimates of p satisfy Equation (13), 
where F,, , F,, and F,, are determined from Equations (14), (15) and (16). 
Consequently Equation (10) must include, as particular cases, the efficient 
estimates given by Stevens. As these estimates are efficient they can 
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also be obtained by minimising the variances of the expression given in 
(10). 
It is of interest to note that, for a given r, 


where the J are the polynomials used by Pimentel Gomes in Equation 
(1). 

Values of yu, , uw. --- appropriate to various values of p in the range 
p = .3 to p = .5 have been determined for the cases n = 5, 6 and 7, 
using the F,,, , F,, and F,, provided by Stevens. The following values 
have been found to lead to estimates of high efficiency (over 90%) over 
the practically useful range of p:— 


n= 5: = 1.75, = 1.50. 

n=6: = 2, = 2.5, ps = 1.75. 

= 7: = 2, = ps = 3, 2 
The actual estimates are | 


= t+ — yo — 
+ 3y2 — yi — 


= + 2ys — 3y2 — 
dy. + 4ys + 2y2 — — Tyo 


Yot Ys + Ys — Y2 — f = 7 19 
Ys + t Ys — — 2Yo (19) 


In the case n = 7 slightly better results can be obtained by using 
Mw, = 2.2, uw. = 3.1, us = 3.3 and uw, = 2.2. These values give almost 
fully efficient estimates when p is about 0.4. The values given above do 
not give efficient estimates for any p but have the advantage of being 
_ whole numbers. 

The estimates (8), (17), (18) and (19) are likely to be most useful 
when values of p are required from a number of similar experiments, 
possibly to check on an assumed value. Thus many problems of deter- 
mining optimum dressings of fertilizer can be solved satisfactorily by 
using the values of k = —log,. p determined from large numbers of 
experiments and suggested by Crowther & Yates (1941). In such cases 
the assumed values of p can be rapidly and easily checked. Full effi- 
ciency is not required for such purposes. 

The proposed estimates can also be used to advantage as preliminary 
estimates of p in Stevens’ efficient method, when the computations 


for n= 5, (17) 


for n = 6 (18) 


and 
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have to be carried out on a desk machine. In most cases only a single 
iteration is required. 

Example 


As an example of the procedure when n = 5 the data given by Stevens 
(1951), p. 261, will be considered. The y are total yields of wheat in 
kgs over five plots of a Latin square: 


y: 44.4 54.6 63.8 65.7 68.9 Mean 59.48 
=: 0 1 2 3 4 
The estimate of p is simply 


, = 1X 68.9 + 3 X 65.7 — 63.8 — 6 X 54.6 
4X 65.7 +3 X 63.8 — 54.6 — 6 X 44.4 


Values of r* are: 


es 1 0.610 0.372 0.227 0.138 Mean 0.4694 


The regression of y on 1” is 


= 0.610. 


This is an estimate of —8. a is estimated by 
59.48 + 28.66 X 0.4694 = 72.93. 


These estimates show close agreement with the efficient estimates which 
Stevens gives as: 


= 7243 45.57, B=28.2445.41, p= 0.597 + 0.140. 


When used as a preliminary estimate in Stevens’ method, r = 0.610 
leads to the following estimates of a, 8 and p:° 


a = 72.38, B = 28.21, p = 0.595. 
Efficiencies of the proposed estimates 


Stevens showed that the variances of the efficient estimates of p are 
given by 


F,,0° 


The variances of the estimates given above can be obtained by applica- 
tion of the usual formula for the variance of a ratio: 


° 
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When applied to the efficient estimates this formula gives results identical 
with those provided by Stevens. 
Thus the efficiency of (4y; + y2 — 5y:)/(4y2 + y: — 5yo) is 


F,,(4p" + p — 5)” 
{42(1 + + 2p} 


The percentage efficiencies of the estimates (8), (17), (18) and (19) 
are set out in Table 1. 

As might be expected the efficiencies fall off for low and high values 
of p. This is, however, of little consequence since for these values 

r'/* is in any case inaccurately determined whether an efficient estimate 
of p is used or not. 

Thus, for n = 7 if s is such that p = 0.1 or 0.7 the variance of r 
is about three times as great as when p = 0.4. In general high and low 
values of r indicate that the intervals between levels of fertilizer are 
too small and too great respectively for the accurate determination of 
curvature. 

The two series of estimates which tend to be fully efficient as p tends 
to 1 or 0 can, however, be expressed simply and may be useful occasion- 
ally. The first series, which is more efficient than the estimates pro- 
posed above for p greater than about 0.6, is given by substituting the 
values wu; — uy, = (¢ + 1)(n — 7 — 1)(n — 2 — 2)/(n — 1)(n — 2) in 
Equation (10). The second series, for which »; = 7 + 1, covers the 
range of p less than about 0.15. 


1/s 


SUMMARY 
Simple methods for fitting the regression equation 
y Bp, where 0<p<1l 


are suggested for cases in which the ordinates x are equally spaced and 
take 4, 5, 6 or 7 values. 

The estimates of p are given by the ratios of two contrasts between 
the values of y. These contrasts are chosen so that the efficiencies of 
the proposed estimates are high over the range of useful values of p. 

Once p has been estimated a and 6 can be determined by a straight- 
forward linear regression. 
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INTERNATIONAL BIOMETRIC SYMPOSIUM ON “THE ROLE 
OF BIOMETRIC TECHNIQUES IN BIOLOGICAL RESEARCH” 


TECHNICAL PROCEEDINGS 


Institute de Educacao Carlos Gomes, Campinas, Brazil, July 4-9, 1955 
Other papers will appear in the next issue 


SOME PROBLEMS OF EXPERIMENTAL DESIGN AND 
TECHNIQUE WITH PERENNIAL CROPS* 


8. C. Pearce 
East Malling Research Station, Kent, England 


Most statistical problems of experimentation with perennial plants 
arise from the longevity of the material, though a few are the result of 
the large size of the plants. Satisfactory classification is not easy, but 
here the problems will be considered in three groups. 

1. Problems in scientific approach arising because few experiments 
are possible in a reasonable period of time, and consequently the best use 
must be made of them. 

2. Problems in trial design and interpretation, mostly of a mathe- 
matical nature, concerning desirable developments ih statistical tech- 
niques. 

3. Problems in the variability and measurement of plants, and the 
relations between various measurements. 


1. Scientific approach 

When experiments each last for several years it is not possible to 
do many of them; but there are sciences, geology and astronomy, for 
example, in which no experiments are possible at all, yet valuable 
advances in knowledge have been made. The difficulty of experimenting 
with long-lived plants does not mean therefore that the acquisition of 
scientific knowledge is impossible; but it does call for a consideration 
of the means available and a careful study of their limitations. 


*Paper presented to the International Bi tric Symposium, Campi Brazil, on 5th’ ‘July 1955. 
A somewhat longer version, similar in substance though not in detail, was kindly prepared by Senhor 
Ediliberto Amaral from an original by the writer, and entitled ‘‘Problemas especificos de delineamento 
e técnica experimental em culturas perenes’’. 
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The obvious alternative is to make good use of systematized obser- 
vations, as is done in the sciences mentioned. An astronomer cannot 
command an eclipse of the sun, but he can forecast when the next will 
take place and study that; in like manner, it might be thought, a research 
worker, unable to carry out a trial on, say, the irrigation of orange trees, 
might examine plantations in which irrigation had taken place and 
compare them with others that had not been so treated. With annuals 
this method has worked very well. It might be expected to work with 
trees also, but it must be noted that many additional factors are in- 
volved. Thus, it is not uncommon to find fields of annuals all of the 
same variety treated in much the same way except for the factor under 
study, the plants being necessarily of the same age; but plantations of 
perennials can differ in so many respects, such as age, rootstock, tree 
form and spacing, that in survey methods it is rarely possible to disen- 
tangle all the factors. 

The same consideration applies to extensive experiments carried out 
on a number of commercial plantations. If similar results are obtained 
at each site, the agreement is impressive, and any recommendation 
based upon it must be of wide application; but if some sites give dis- 
cordant results it is rarely possible to identify with certainty the factor 
responsible. 

Nevertheless, much knowledge of perennial plants has been obtained. 
Because of the difficulties of experimenting, however, much of it has 
come by other channels, and it is necessary to recognize this and to 
consider if these other channels can be improved. One valuable source 
of knowledge, for example, is the advisory officer. When visiting 
growers he receives as.well as gives information, and in this way he 
accumulates evidence until it is clear to him either that a certain new 
method is an improvement or that opinion about it is so divided that 
it cannot represent an important advance. Such an approach can 
readily be systematized statistically; a group of growers selected at 
random from those with experience of the matter might be approached, 
and if their replies were inconclusive a further group might be asked, 
and so on until a result was obtained. Since much knowledge, especially 
about new varieties, does in fact come by this channel, it would be a 
good thing if it could be obtained more systematically and with proper 
statistical safeguards. 

Basically, however, the problem must be to improve experimenta- 
tion, and the best way of doing this is to have a clear purpose behind 
each trial. In any large set of treatments it is apparent that many 
comparisons can be made between them or between groups of them. 
With only six treatments, for example, there are 301, and by studying 
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enough comparisons something significant can usually be found. In 
practice this does not matter, because no one does test separately all 
possible comparisons but only those that have an intelligible meaning; 
but to select these requires some knowledge outside that given by the 
trial itself. 

Reflection will show that this outside knowledge is used in three ways: 

‘a. by’ previous explanation. An experimenter by studying existing 
knowledge perceives certain laws or regularities, and wishes to test if 
these are coincidental or genuine by carrying out an experiment, de- 
signed to give a specified result if he has been thinking on the right lines, 
but not otherwise. 

b. by subsequent explanation. An experiment has been performed 
the result of which was not foreseen, but later consideration suggests 
an explanation. 

c. by further examination. A result has been obtained for which no 
explanation can be advanced. In view of its importance, if true, further 
work is done, and the same result appears again. 

Of these three, the last is mere empiricism, the gathering together 
of facts, not the understanding of the permanent elements that lie 
behind facts, which is the domain of science. With annual plants, 
where experiments are readily repeated, it provides a useful practical 
approach. With long-lived plants, on the contrary, where it is more 
difficult to accumulate facts by sheer weight of evidence, it is necessary 
always to seek an understanding of what is being observed. But this 
understanding should be sought as early as possible. The second of 
the ways mentioned above, obtaining a result and thinkjng up an explan- 
ation afterwards, is rarely convincing; if the explanation is complex it 
is judged to be over-ingenious and if simple to be belated. At the best 
it is a chancy method, and depends very much on the integrity of the 
experimenter. The most powerful experiments are those in the first 
class, for they do not only establish a fact but justify a system of thought. 
When Galileo dropped two balls from the Leaning Tower of Pisa 
he did not merely show that gravitational acceleration is independent of 
mass but he shattered the whole structure of Aristotelian physics. 
Where experiments are necessarily few and costly it is powerful ones 
like these that must be sought. 

Also, if the end of experimentation with perennial plants must 
always be comprehension, it is not enough merely to measure crop or 
whatever is immediately under study. It does not suffice to assert 
that in a particular trial Treatment A gave more crop than B without 
‘knowing why, for there will rarely be enough empirical evidence to 
establish the fact as a general rule, and its broad acceptance must 
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derive from the mechanism of the effect being inherently reasonable 
and of a sort likely to operate in a wide range of conditions. Conse- 
quently, the good experimenter will measure not only the weight of 
fruit, but also growth in its various aspects, disease, leaf colour, amount 
of blossom, and perhaps much else, so that his crop figures do not stand 
in isolation, but are seen as the culmination of a biological process, 
each part of which has taken place beneath a comprehending scrutiny. 


2. Trial design 

Problems in trial design that arise with perennial plants are not at 
the moment especially numerous, possibly because of the relatively 
large amount of work that has already been done. 

First, there is the need to develop row and column designs so that 
best use can be made of limited spaces. Perennial plants are mostly 
large, so if an area will take only seven rows of five trees there is room 
for no more than 35 plots; it is no use complaining that with six treat- 
ments 36 plots would be much more convenient. Further, if efficient 
use is to be made of the area, it is desirable to adopt a design that will 
allow for positional effects by the use of rows and columns rather than 
by blocks, for such designs will allow for the differences between outside 
and inside trees and so permit both to be used, subject to certain obvious 
restrictions. It happens that the evolution of such designs is formally 
part of a larger problem that at first sight is quite different, namely, 
that of changing the treatments on long-lived plant material. Whenever 
an experiment is concluded while the plants are still in their prime, as 
often happens, the possibility of using the same trees for further experi- 
ments must be examined. If the first trial is in randomized blocks, the 
problem is the same as that already mentioned, namely, that of adding 
a third classification to plots already classified in two mutually ortho- 
gonal ways (rows and columns in the one instance, blocks and original 
treatments in the other), and of doing so in such a manner that the 
new classes are compared to greatest advantage. Of course, the first 
trial is not necessarily in randomized blocks, so the problem of changing 
treatments is the wider one, and in its mathematical aspects includes 
that of using restricted areas economically. The complete solution of 
this problem must await the compilation of a comprehensive catalogue 
of balanced designs, but already considerable progress has been made. 
Thus, in the problem mentioned, i.e., applying six treatments to 35 
trees arranged seven by five, the modern experimenter has the resource 
of using a Latin square with a column added and a row omitted. 

The large size of most perennial plants suggests the use of parts of 
an organism, e.g., the branches of a tree, as plots of an experiment, the 
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word “organism” underlining the possibility of treatments so applied 
affecting other plots of the same block. Thus, in an experiment upon 
pollination some blossoms of a tree may be pollinated with Variety A 
and others with B; but it may be that the good set of the first will inhibit 
setting elsewhere, including those blossoms pollinated with B. A little 
reflection will show that given complete block designs there is no way 
of distinguishing the local effect, 7, of a treatment upon the plots to 
which it is applied from the remote effect, 5, upon the plots elsewhere 
in the organism. The plots to receive it show, for example, ten per 
cent better set than those that have been treated otherwise, but is this 
because they have been improved, or the rest have been worsened, or 
something of each? The difference equals (y — 5) and there is no 
means of determining the individual values of the two parameters. If, 
however, incomplete block designs are used, the treatment effects shown 
by intra-block comparisons equal (y — 64) as before, but those from 
inter-block comparisons equal [fy + (k — 1)é], where k is the number 
of plots to a block. This technique is already being used in an exploratory 
way at East Malling, where much of the mathematics has been worked 
out. Thus, experiments have been carried out on the effects of different 
kinds of pollen on fruit set and on the effects of different degrees of 
defoliation on fruit size. There is advantage both in the use of plots 
smaller than a whole tree and in the avoidance of doubtful assumptions 
about relationships within an organism. Further, isolation of the remote 
effects sometimes sheds light on the mechanism of the treatments under 
study. 

The last problem in this group that will be considered concerns the 
combining of several years’ results. It is dangerous to assume that the 
total crop over a period is all that matters, and equally dangerous to 
concentrate on single years in isolation; the pattern of bearing, whether 
the cropping tends to increase, or to decrease, or to fluctuate, or to be 
irregular or regular, may be important. To a large extent these charac- 
teristics are determined by the varieties themselves and the age of the 
trees, but they also may be affected by the treatments and it is most 
important to know to what extent. Probably the best way of dealing 
with this problem is one associated with Brazil, namely, that suggested 
by Professor Stevens of Séo Paulo, in which various linear functions of 
annual yields are examined successively, each function being chosen to 
measure some characteristic of the pattern of cropping, such as the 
biennial swing, or the shape of the curve after this has been allowed for. 
The method can be used equally well for growth, incidence of disease, 
or anything else. One proviso may be mentioned here. The fitting 
of curves to past crops is a useful descriptive device; their extrapolation 
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to give information about future ones is quite another matter. The 
only satisfactory way of finding out about crop during twenty years 
is to grow the trees for that period. 

Valuable though this method is—it has incidentally the advantage 
of not being vitiated by the variability differing from year to year—it 
does not go far enough. For one thing, cyclic effects are not always 
biennial, and it has been suggested that for some crops, such as coffee, 
the period is between two and three years, and with forest trees much 
longer. Also, a biennial rhythm can be thrown out of phase by meteoro- 
logical factors. For another thing, the general shape of the curve after 
smoothing out cyclic effects is not easily expressed in any way that 
makes good biological sense. The general approach is no doubt sound, 
namely, the evolution of tests to supplement an analysis of variance on 
total crop borne over the whole period; but these tests, when obtained, 
may well prove more subtle than examination of linear functions. 
Here is a problem that needs much more thought than it has received 


' up to the present. 


3. Variability and measurement 


Although the second group of problems may be said to be few and 
simple, those of the last group are numerous and complex, and also 
have to be considered anew for each species. 

One of the first importance is the determination of sources of error, 
because this may largely affect the design. Variation can come from 
two main sources, the trees themselves or positional effects, and when- 
ever a trial is proposed it is necessary to be clear if one or both of these 
sources needs to be considered, otherwise a lot of trouble can be taken 
controlling variation that is negligible anyway. 

In particular, if most variation arises in the plants themselves there 
will be a great need to find characteristics apparent at planting time 
that can be used to give worth-while adjustments by the method of 
covariance. There would not, however, be much gained by attempts to 
control positional effects, and any complications of design introduced 
for this purpose would probably on the whole do more harm than good. 
On the other hand, if the variation is mostly positional, devices such as 
confounding and incomplete blocks may be valuable, and it would also 
be worth-while to go to the trouble of using many small plots rather 
than a few large ones. 

When trees come from the nursery for planting it is advisable to 
put them into their permanent places at random. Where this is done 
the initial variation is associated with the trees, and has no positional 
elements in it. Then, when the trees begin to grow, further variation 


eve, 
) 
e 
st 
£ 
of 
0 
rr. 
e, 


336 BIOMETRICS, SEPTEMBER 1956 


is introduced on account of differences in growth rate induced either by 
position, or by genetical and other factors inherent in the trees. In 
fact, some positional variation may well arise, but it may be negligible 
or it may dominate according to circumstances. 

The ideal is, no doubt, to use plant material of the greatest uni- 
formity—and it is impossible to stress too much the need for good, 
standardized plants—and also to use the smallest blocks possible so as 
to minimize positional variation. However, a certain sense of proportion 
is called for. Errors combine in such a way that reduction of the smaller 
component is rarely worth bothering about; it is always the larger 
that must be tackled. But which, for any species in any locality, is the 
larger? So often no one has ever thought to enquire, but it is a basic 
problem. Of course, initially the variation always comes from the 
trees; the question should be put “Will the positional effects come to 
matter?’’. 

Other biometric problems are those of measurement and the rela- 
tionships between measurements. When discussing questions of scien- 
tific approach to investigations with perennial plants it was suggested 
that a comprehending scrutiny should be maintained; it is all very well 
to say that, but how is it to be done? For example, the activities of a 
woody plant depend upon its ability to form a structure from which 
leaves can be displayed to the sun’s rays, but has anyone ever found 
out how to measure this extension growth rapidly but effectively? On 
the other hand, the trunk circumference can be measured easily enough, 
but does anyone know what it means? In general, of course, the thicker 
the trunk the larger the tree; but with apples, and probably with other 
species also, severe pruning leads to growth which at first finds no reflec- 
tion in the trunk. Also, in apples, the relationship of tree weight to 
trunk circumference depends upon the rootstock; so evidently the 
measurement must be used with some circumspection in trials of 
pruning methods and of rootstocks. Again, growth is apparently 
inhibited by cropping, which is related to blossoming, but is it the crop 
or the presence of fruit buds that inhibits the growth? There is some 
evidence that it is the blossoms that check the growth; but comprehen- 
sion of what is happening in a tree is impossible until fairly clear answers 
have been found to this and many similar questions. Admittedly 
these are physiological problems also, but the statistician with his 
methods of regression analysis and variance components has no excuse 
for idleness either. It is not too much to say that inadequacy of measur- 
ing techniques and the ignorance of relationships between measure- 
ments are the chief weaknesses in present-day experimentation with 
perennial plants, and ones that present a continual challenge to the 
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statistician. It has been suggested in the first section that an experiment 
must test a way of thought, that it must be a seeking of comprehension 
of the factors behind phenomena—that empiricism is not enough; but 
this is futile if no means are at hand for exploring what goes on inside 
a tree without killing it to find out. The need all the time is to wnder- 
stand, and this we often cannot do because we cannot measure. 

In putting forward these problems for consideration no attempt has 
been made at comprehensiveness, although those advanced are perhaps 
fairly representative of the range and variety of problems to be en- 
countered. Their enumeration will have served a useful purpose if it 
does no more than show that experimentation with perennials is a 
branch of statistics with a number of problems of its own, some of which 
have received far too little attention. 
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QUERIES 


GeorcE W. SNEpEcoR, Editor 


QUERY: I would appreciate your advice regarding the statistical 
122 analysis of the following problem. The problem concerns the 

length of gestation period in cattle and the data are classified 
according to three factors (1) sex of the calf born (2) the pregnancy 
order of the dam (3) month of calving. This is then a2 X R X C Table 
where R = 7 and C = 12 in the present data. The data are coded by 
subtracting 260 days from each gestation period (the range was 264 
to 296 days). The number of observations in the cells of the table 
varied from 0 to 19, there being 15 cells with no observation. Values 
were calculated (to the nearest whole number) for these 15 cells from 
the means for sex, pregnancy number and month of calving on the 
assumption that there was no interaction between the various factors. 

The information I wish from the data is as follows: 


. The difference in gestation length because of the sex of the calf. 

. Differences in gestation length due to pregnancy order. 

. Seasonal differences in gestation length. 

. The presence or absence of interactions between these factors. 

. Estimates of the proportions of total variance attributable to these 
factors. 


ork Whe 


Had no cells been missing in your 2 X 7 X 12 table I 
ANSWER: would have recommended an unweighted analysis of 

the cell means except for intracell variation although 
the number of observations per cell varies considerably. Missing cells, 
however, and particularly as many as 15 missing cells, necessitate a 
least squares analysis in order to obtain the itemized information, if 
the full classification of each factor is to be maintained. But the in- 
clusion of interactions in the model makes this metiicd too laborious 
and impractical, if not impossible. 

A feasible alternative, when one or more of the classifications of a 
factor can be logically grouped together so that no cells are missing, 
is to form a new 2 X R’ X C’ (R’ < R and/or C’ < C) table with no 
missing cells. This procedure should be applicable to your data since the 
missing cells most likely come from having few observations on dams 
with high pregnancy orders (5, 6 and 7). It seems probable, with only 
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15 cells missing out of 168 as you now have the data classified, that 
grouping pregnancy orders 6 and 7 together or at most 5, 6 and 7 
together would leave no missing cells. Seasonal calving would lead to 
few observations in certain months and this could be overcome by 
grouping the months into seasons. However seasonal calving is not 
generally regulated to the extent that it is expected to be a prominent 
factor in disproportionality. 

Having grouped the classifications so that no cells are missing the 
usual analysis of variance of the cell means can be performed. The 
expectations of the mean squares are those for a single observation in 
each cell except that the coefficient of the within cell variance is the 
reciprocal of the harmonic mean number of observations in the cells. 
The within cell variance is computed separately and pooled over cells. 

The information you wish through item 4 is obtained from the above 
analysis of variance or by direct comparisons of unweighted means. 
For item 5, to estimate the proportions of total variance attributable to 
these factors, it must be remembered that the model corresponds to a 
fixed one in that the classifications cover the full range of effects for sex 
and month of calving and for all practical purposes for pregnancy order. 
However, I can see no practical harm in assuming a completely random 
model and to apportion the total variance to the various sources by this 
method. 

The grouping of adjacent classes where an order is defined by the 
classification such as in pregnancy order reduces the variances stemming 
from this source. However, when the classes grouped are not very 
different, and in the case of pregnancy order, the orders 5, 6 and 7 are 
usually about the same for most characters, the reduction in variance 
should be slight, if any. 


C. CLarK CocKERHAM 


CORRECTION TO “FRACTIONAL REPLICATION FOR 
MIXED SERIES” 


Mitton Morrison 
Experimental Towing Tank, Stevens Institute of Technology 


Dr. K. A. Brownlee has pointed out that there is a slip on page 11 
of my paper (this journal, Vol. 12, 1956, pp. 1-19). The degrees of 
freedom for E, AE, BE, CE and DE should each be 3, and not 2. Those 
for the residual should be correspondingly reduced from 11 to 6. 


tts 
4 
! 
3 
G ft 
Sy 
belt 
AG 
wit 
} 
> 


NEWS AND ANNOUNCEMENTS 


Members are invited to transmit to their National or Regional Secretary 
(if members at large, to the General Secretary) news of appointments, 
distinctions or retirements and announcements of professional interest. 


DR. JOHN WISHART 


Dr. John Wishart, Reader in Statistics at Cambridge University, 
drowned by misadventure while swimming from El Revolcadero beach, 
Port of Acapulco, Mexico on July 14th, 1956. Dr. Wishart was dis- 
tinguished for his theoretical contributions, especially in respect of 
the distribution of moment and product moment statistics, and for his 
applied work in agricultural experimentation. He was a charter member 
of the Biometric Society. He had been in Mexico since January, 1956 
in connection with the Centro Interamericano de Capacitation sobre 
el Uso de Metodos Estadisticos en la Experimentation Agricola, spon- 
sored by the Government of Mexico, and by FAO. 


Professor Ralph A. Bradley was the recipient of the Brumbaugh 
Award of the American Society for Quality Control, given annually 
to the author of the article published in Industrial Quality Control, 
judged to have made the greatest contribution to the development of 
industrial applications of quality control. 


Under an agreement into which the Johns Hopkins University 
entered some years ago, Professor W. G. Cochran was visiting professor 
at the University of the Philippines during the months of June and 
July 1956. 


The eminent distinction of the Copley Medal of the Royal Society 
of London was awarded to Professor Sir Ronald A. Fisher, F.R.S., 
‘fn recognition of his distinguished contributions to developing the 
theory and application of statistics for making quantitative a vast 
field of biology’. 


Dr. C. H. Goulden, previously awarded the Gold Medal of the 
Professional Institute of the Public Service of Canada for his scientific 
contributions, is now Director, Experimental Farms Service, Canada 
Department of Agriculture. 


Professor George W. Snedecor’s many outstanding contributions to 
statistics were recognized when a honorary Doctorate of Laws was 
conferred on him by the University of North Carolina, May 24, 1956. 
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THE BIOMETRIC SOCIETY 
International Affairs 


The Society has been asked to organize three sessions for the meeting 
of the International Statistical Institute to be held in Stockholm from 
August 8 to 15, 1957. Plans are also being made for the 4th Inter- 
national Biometric Conference to be held in Ottawa during 1958, 
probably from August 28 to September 1. Details of these meetings 
will be sent to members as soon as available. 

Biometric Seminars are being held in Linz from September 24 to 
October 3, 1956, and in Milan from October 8 to October 20. 


Switzerland 


The Swiss section of the Society held its annual meeting on July 
14, 1956 in Zurich. The programme included communications by H. L. 
LeRoy: Die Inzuchttheorie von R. A. Fisher; P. Geier: Les rapports 
numériques entre un Lépidoptére nuisible aux cultures fruitiares et 
certains de ses antagonistes. 

At this meeting, Dr. H. L. LeRoy was elected National Secretary 
in place of Professor A. Linder. 


ENAR 


A joint meeting with the American Institute of Biological Sciences 
was held at Storrs, Conn. on August 28-29. The programme included 
the following papers:— 

Joint session with the Ecological Society of America—L. C. Cole: 
Some examples of the ecologist’s problems; C. I. Bliss: Some over- 
dispersed distributions; P. J. Clark: Nearest Nees methods in the 
analysis of spatial pattern. 

Session of contributed papers—C. White: Distribution of erythro- 
cytes on the floor of the hemocytometer chamber; I. E. Satterthwaite: 
Random assignment experimental designs. 

Joint sessions with the American Society for Horticultural Science— 
R. L. Wine: On the comparison of multiple test procedures; C. W. Dunett: 
Multiple decision procedures. 

A joint meeting with the American Statistical Association was held 
in Detroit, Mich. from September 7-10. The programme included:— 

*Covariance Analysis—D. B. deLury: Elements of covariance; 
W. T. Federer: Covariance analysis with unequa! subclass numbers; 
M. Zelen: The analysis of covariance for incomplete block designs; 
H. O. Hartley; Group comparisons and analysis of variance and covari- 
ance in cluster sampling. 
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Animal Feeding Experiments—G. C. Ashton: Control of error 
variance in swine feeding experiments. 

Screening problems in industrial biology—C. W. Dunnett: The 
Pharmaceutical Industry; M. B. Carroll: The Food Industry; H. Smith, 
jr: Biochemical Research. 

*Applications of stochastic processes—D. W. Alling: The after- 
history of pulmonary tuberculosis; A. F. Bartholomay: The kinetics 
of enzyme action. 

Biological standardization—I. A. DeArmon: Control of precision 
in the plate count assay; B. W. Brown: More sensitive vaccine screening 
procedures; J. Ipsen: Assays of immunizing agents; M. C. Sheps: 
Between-assay error in biological assays. 

The sessions marked * were co-sponsored by the Institute of Mathe- 
matical Statistics. 


German Region 


A Biometric Colloquium at Hanover is planned for September 
10th, with the following programme—A. Lein: Anwendungsmdglich- 
keiten der Biometrie in der Landwirtschaft; M. Wermke: Die Ver- 
wendung von Lochkartenmaschinen zur Auswertung von Feldversuchen; 
W. Seyffert: Versuchsauswertung nach Messung mehrerer Merkmale; 
H. Rundfelt: Zur Rationalisierung der Pflanzenziichtung. 


WNAR 


The Annual Meeting was held at Seattle on August 22-24. Included 
in the programme was an invited address by Prof. N. Rashevsky, 
entitled ““Models and General Mathematical Principles in Biology and 
Sociology”’. . 

Other papers given at the meeting included:—L. Moses: Some 
non-parametric techniques; W. Becker: The log transformation of 
growth data; D. Jenden: Reaction rates in geometrically constrained 
enzyme systems; H. Hotelling: New light on the multiple correlation 
coefficient; P. Horst: Optimal estimates of multiple criteria with restric- 
tions on the covariance matrix of estimated criteria; V. Miller: Procedural 
consideration in forecasting population; F. Johnson: The role of sta- 
tistical methods in forest and range research; T. Orr: Some notes on 
the growth of partially cut stands; J. H. G. Smith and J. W. Ker: 
Some distributions encountered in sampling of -forest stands; B. Wagg: 
The relationship between the growth of the Spruce budworm and 
certain environmental factors; M. Ahmed: A stochastic model for the 
tunneling and re-tunneling of flour beetles; E. R. Rich: A stochastic 
model for the number of beetles on the surface of flour. 
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BIOMETRICS—DIRECT INCOME AND EXPENSE—1955 


(Final Report, March 29, 1956) 
Income 
(1) Reserve (1951 thru 1954) sale 3,1 849.00 
Reserve (1951 thru 1954) sale 7,1 756.00 


Balance from 1954 books 6387 .34 
7992.34 
(2) Member subscriptions 
ASA 443 at 4.00 1772.00 
Bio. Soc. 573 at 4.00 
680 at 2.75 
6 at 25.00 4312.00 
6084 .00 
(3) Non-members subscriptions (690) 4835.75 
(4) Sale of back issues 
Vol. 3, No. 1 106 at 1.50 159.50 
Vol. 7, No.1 23 at 2.00 46.00 
Others 2038 . 47 
2243 .97 
(5) Sale of reprints 998 . 67 
(6) IUBS (Proceedings Brazil 
Symposium) 500.00 
22654 .73 
Expenses 
(1) Overpayment and cancelation 
of subscriptions 25.00 
(2) Joseph Ruzicka, bookbinding 23.25 
(3) Express charges 598.62 
(4) Inst. of Stat. editorial manage- 
ment 1000.00 
(5) Supplies 239.69 
(6) Postage 300.00 
(7) J. W. Hopkins 500.00 
(8) Wm. Byrd Press, reproduction 
of Vol. 9, No. 1 953 .96 
(9) Wm. Byrd Press Production | 
of Vol. 11, No. 1 1799.46 
Vol. 11, No. 2 2066 . 33 
Vol. 11, No. 3 2234.70 
Vol. 11, No. 4 2732.13 
8832.62 
Offprints of 
Vol. 10, No. 4 337.31 
Vol. 11, No. 1 286 .07 
Vol. 11, No. 2 231.73 
Vol. 11, No. 3 365.09 
Vol. 11, No. 4 358.94 
1579.14 
10411.76 
14052. 28 
Balance 
8602.45 
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